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PREFACE 



This book is the result of an attempt to answer a number 
of questions which frequently confront engineers. So far as 
the author is aware no other book in English covers the same 
groimd in an elementary manner. 

It is thought that the method of determining the constants 
in formulas by the use of the straight line alone leaves little to 
be desired from the point of view of simplicity. The approxi- 
mation by this method is close enough for most problems arising 
in engineering work. Even when the Method of Least Squares 
must be employed the process gives a convenient way of obtain- 
ing approximate values. 

For valuable suggestions and criticisms the author here 

expresses his thanks to Professors Alexander Ziwet and Horace 

W. King. 

T. R. R. 

University of MicmoAN, 1917. 
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INTRODUCTION 

In the results of most experiments of a quantitative nature, 
two variables occur, such as the relation between the pressure 
and the volume of a certain quantity of gas, or the relation 
between the elongation of a wire and the force producing it. 
On plotting the sets of corresponding values it is foimd, if they 
really depend on each other, that the points so located lie 
approximately on a smooth curve. 

In obtaining a mathematical expression which shall represent 
the relation between the variables so plotted there may be two 
distinct objects in view, one being to determine the physical 
law underlying the observed quantities, the other to obtain 
a simple formula, which may or may not have a physical basis, 
and by which an approximate value of one variable may be 
computed from a given value of the other variable. 

In the first case correctness of form is a necessary considera- 
tion. In the second correctness of form is generally considered 
subordinate to simplicity and convenience. It is with the 
latter of these (Empirical Formulas) that this volume is mostly 
concerned. 

The problem of determining the equation to be used is really 
an indeterminate one; for it is clear that having given a set of 
corresponding values of two variables a number of equations 
can be found which will represent their relation approximately. 

Let the coordinates of the points in Fig. i represent different 
sets of corresponding values of two observed quantities, x and y. 
If the points be joined by segments of straight Unes the broken 

9 
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line thus formed will represent to the eye, roughly, the relation 
between the quantities. 

It is reasonable to suppose, however, that the irregular dis- 
tribution of the points is due to errors in the observations, and 
that a smooth curve drawn to conform approximately to the 
distribution of the points will more nearly represent the true 
relation between the variables. But here we are immediately 
confronted with a difficulty. Which curve shall we select? a 
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Fig. I. 



or J? or one of a number of other curves which might be drawn 
to conform quite closely to the distribution of the points? 

In determining the form of curve to be used reliance must 
be largely placed upon intuition and upon knowledge of the 
experiments performed. 

The problem of determining a simple equation which will 
represent as nearly as possible the curve selected is by far the 
more difficult one. 

Ordinarily the equation to be used will be derived from a 
consideration of the data without the intermediate step of 
drawing the curve. 

Unfortunately, there is no general method which will give 
the best form of equation to be used. There are, however, a 
number of quite simple tests which may be appUed to a set of 
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data, and which will enable us to make a fairly good choice of 
equation. 

The first five chapters deal with the application of these 
tests and the evaluation of the constants entering into the 
equations. Chapter VI is devoted to the evaluation of the 
constants in empirical formulas by the Method of Least Squares. 
In Chapter VII formulas for interpolation are developed and 
their applications briefly treated. Chapter VIII is devoted to 
approximate formulas for areas, volumes, centroids, moments of 
inertia, and a number of examples are given to illustrate their 
application. 

Figs. I to XX at the end of the book show a few of the forms 
of curves represented by the different formulas. 

A few definitions may be added. 

Arithmetical Series. A series of numbers each of which, 
after the first, is derived from the preceding by the addition 
of a constant number is called an arithmetical series. The con- 
stant number is called the common difference 

6, 6.3, 6.6, 6.9, 7.2, 7-5 •• • 

and 

18.0, 15.8, 13.6, 11.4, 9.2 .. . 

are arithmetical .series. In the first the common difference is 
.3, and in the second the common difference is —2.2. 

Geometrical Series. A series of numbers each term of which, 
after the first, is derived by multiplying the preceding by some 
constant multiplier is called a geometrical series. The constant, 
multiplier is called the ratio. 

1.3, 2.6, 5.2, 10.4, 20.8, 41.6 . . . 

and 

100, 20, 4, .8, .16, .032 . . . 

are geometrical series. In the first the ratio is 2, and in the 
second it is .2. 

Differences are frequently employed and their meaning can 
best be brought out by an example. 
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X 


y 


Ay 


A=y 


A'y 


A^ 


I 


I0.2 










2 


II. I 


0.9 
I.I 


0. 2 


0.0 




3 


12.2 


1-3 

2.7 
1.8 


0.2 


1.2 


1.2 


4 
S 
6 


I3S 
l6.2 

i8.o 


1-4 
-0.9 
-0.8 


-2-3 

O.I 


-3-5 
2.4 
2.7 






I.O 




2.8 




7 


19.0 




2.0 






8 


22.0 


3-0 









In the table corresponding values of x and y are given in 
the first two columns. In the third column are given the values 
of the first differences. These are designated by Ay. The first 
value in the third column is obtained by subtracting the first 
value of y. from the second value. The column of second differ- 
ences, designated by l^y, is obtained from the values of Ay in 
the same way that the column of first differences were obtained 
from the values of y. The method of obtaining the higher differ- 
ences is evident. 



CHAPTER I 

I. y = a+bx+cx^+dx'^+ . . . +q3f 

Values of x form an arithmetical series and A"y constant. 

In a tensile test of a mild steel bar, the following observa- 
tions were made (Low's Applied Mechanics, p. i88): Diameter 
of bar, unloaded, 0.748 inch, ]F=load in tons, x = elongation in 
inches,'On a length of 8 inches. 



w 


I 


2 


3 


4 


5 


6 


X 

Ax 


0.0014 
0.0013 


0.0027 
0.0013 


0.0040 
0.0015 


0.005s 
0.0013 


. 0068 
0.0014 


0.0082 



Plotting W and x, Fig. 2, it is observed that the points lie 
very nearly on a straight hne.* Indeed, the fit is so good that 
it may be almost concluded that there exists a hnear relation, 
between W and x. From the figure it is found that the slope 
of the line is 0.00137 ^.nd that it passes through the origin. The 
relation between W and x is therefore expressed by the equation 

x=o.ooi^'jW. 

The observed values of x and the values computed by the 
above formula are given in the table below. 



XV 


Observed x 


Computed x 


1 


0.0014 


0.00137 


2 


0.0027 


0.00274 


3 
4 
S 
6 


0.0040 
0.0055 
0.0068 
0.0082 


0.0041 I 
0.00548 
0.00685 
0.00822 



* By the use of a fine thread the position of the line can be deter- 
mined quite readily. 
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The agreement between the observed and the computed 
values is seen to be quite good. It is to be noted, however, that 
the formula can not be used for computing values of x outside 
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Fig. 2. 



the elastic limit. In the experiment 6 tons was the load at the 
elastic limit. 

It is not necessary to plot the points to determine whether 
they lie approximately on a straight line or not. Consider the 
general equation of the straight line 

y=mx-\-k. 

Starting from any value of x, give to x an increment. A*, and 
y will have a corresponding increment, Ay. 

y-\-Ly = m{x-\-l\o^-\-k\ 
y=mx+k; 
Ay — mAx. 
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From this it is seen that, in the case of a straight line, if the 
increment of one of the variables is constant, the increment of 
the other will also be constant. 

From the table it is observed that the successive values of 
W differ by unity, and that the difference between the successive 
values of x is very nearly constant. Hence the relation between 
the variables is expressed approximately by 

x=mW+k, 

where m and k have the values determined graphically from 
the figure. 

By the nature of the work it is readily seen that the graphical 
determination of the constants wiU be only approximate under 
the most favorable conditions, and should be employed only 
when the degree of approximation required will warrant it. 
Satisfactory results can be obtained only by exercising great 
care. Carelessness in a few details will often render the results 
useless. Understanding how a graphical process is to be carried 
out is essential to good work; but not less important is the 
practice in applying that knowledge. 

In experimental results involving two variables the values of 
the independent variable are generally given in an arithmetical 
series. Indeed, it is seldom that results in any other form 
occur. It will be seen, however, that in many cases where 
the values of the independent variable are given in an arithmeti- 
cal series it will be convenient to select these values in a geomet- 
rical series. 

As a special case consider the equation 

y = 2—^x+x^. 

If an increment be assigned to x, y will have a corresponding 
increment. The values of x and y are represented in the table 
below. Ay stands for the number obtained by subtracting any 
value of y from the succeeding value. A^y stands for the num- 



16 



EMPIRICAL FORMULAS 



ber obtained by subtracting any value of Ay from its succeeding 
value. The values of x have the common difference 0.5. 



X 


°-S 


I.O 


i.S 


2.0 


2-5 


30 


3-S 


4.0 


y 


0.7S 


0.00 


-0.2s 


0.00 


0.7S 


2.00 


3-7S 


6.00 


Ay 


-0.7S 


-0.2s 


0.25 


0.7S 


I- 25 


I-7S 


2.25 




A2jl 


0.50 


0.50 


0.50 


0.50 


0.50 


0.50 







The values A^y, which we call the second differences, are 
constant. 

These differences could equally well have been computed 
as follows: 

y+Ay = 2 — ^(x+Ax) + {x+Axy, 
Ay= —^{Ax) + 2x{Ax) + (Axy, 
Ay+A''y= -3(Aa;) + 2(a;+Aa;)(Aa;) + (Aa;)2, 
A^y= 2(Aa;)2 =0.5 since Aa; = 0.5. 

From this it is seen that whatever the value of Ax (in 
y = 2—^x+x^) the second differences of the values of y are 
constant. 

Consider now the general case where the «th differences are 
constant. For convenience the values of y and the successive 
differences will be arranged in columns. The notation used is 
self-explanatory. 



yi 


Ayi 










y2 




A^yi 










Ay2 




A3yi 




etc., 


ys 




A^ya 




A*yi 






Ays 




A^y2 




etc., 


y4 




A^ys 




A^ya 






Ayi 




ASya 




etc., 


ys 


Ays 


A2y4 








ye 




. . . 
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From the above it is clear that 
y2=yi+^yi, 

= 3'i + 2Ayi+A2yi. 

y4=y3+A3'3, 

=yi+2A3'i+A2yi+A(yi+2A3;i+A2yi), 
=3»i +3Ayi +3A23»i + A^yi. 

yb=yi+^yi, 

= y\ H-SAyi +3A^>'i + A^yi + ACyi H-sAyi +3A2yi + A^yi) 
=yi+4A3'i +6A2yi +4A3yi +A4yi. 

In the above equations the coefficients follow the law of 
the binomial theorem. Assuming that the law holds for yx 
it wiU be proved that it holds for yt+\- 

By hypothesis 

y.=yi+(^-i)Ayi+ ^^~'^^^~^^ A2yi 

|2 
,"(;&- l)('^-2)(/fe-3),„ , , , - 

+.^ iv_^ — '^ ^A^yi+etc (i) 

M 

If this equation is true, then 

y.+i =yi + (^ - i)Ayi + ^^~y^~'^A ^yi 

-f- i il-. — ^^ ■^A^yi + etc. 

13 



+A[yi + (/fe-i)Ayi + ^^ '^J^ '^ A^yi 



>-^)^^-^)(^-^A3,.+etc.] 



, , . , ^(^ — i) .9 , ^(^ — 1)(^ — 2)^, 
= yi+M>'i+-ih iA2y+-i pi^^ '-A^y 



K^-i)(^-2)(^-3) 



A*yi+etc. 
4 
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This is the same law as expressed in the former equation, and 
therefore, if the law holds for y^, it must also hold for jt+i. 
But we have shown that it holds for y^,, and therefore, it must 
hold for ys. 

Since it holds for yt, it will hold for ye. By this process it 
is proved that the law holds in general. 

If now the first differences are constant the second and 
higher differences will be zero, and from (i) 

>'* = 3'i + (A!-i)Ayi. 

If the second differences are constant the third and higher 
differences will be zero, and it follows from (i) that 

|2 

In general, then, if the «th differences are constant 

, /L N. , ('fe-i)('fe-2)^2 , (^-I)(^-2)(*-3)A3■ 
yi = yi + (^ - i) Ayi +-^^ jy ^A^yi ■\-- ■^-^ — ^A3y 1 



The law requires that the values of x form an arithmetical series, 
and hence 

xt.=^x\-\-Qi — -C)Lx\ 
from which follows 

*=^+i (3) 

Aa; 

Substituting this value of k in equation (2) it is found that 
the right-hand member becomes a rational integral function of 
a;,: of the wth degree. Equation (2) takes the form 

yt = a^-hXk-\-cx}^-dx?+ . . . +qx". 

Since % and yt are any two corresponding values of x and y 
the subscripts may be dropped and there results the following 
law: 
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If two ■variables, x and y, are so related that when values of x 
are taken in an arithmetical series the «th differences of the cor- 
responding values of y are constant, the law connecting the variables 
is expressed by the equation 

I y=a+bx+cx^+dx^+ . . . +ga;". 

The nth. differences of the values of y obtained from observa- 
tions are seldom if ever constant. If, however, the wth differ- 
ences approximate to a constant it may be concluded that the 
relation between the variables is fairly well represented by I. 

As an illustration consider the data given on page 131 of 
Merriman's Method of Least Squares. The table gives the 
velocities of water in the Mississippi River at different depths 
for the point of observation chosen, the total depth being taken 
as unity. 



Ah 



Ah 



At surface. 
0.1 depth. 
0.2 

0-3 
0.4 

o-S 
0.6 
0.7 
0.8 
0.9 



3-1950 
3.2299 

3-2532 
3.2611 
3- 2516 
3.2282 
3.1807 
3.1266 
3-0594 
2-9759 



+349 
+ 233 
+ 79 
- 95 
-234 
-475 
-541 
—672 

-835 



-116 

-154 
-174 

-139 
-241 
- 66 

-131 
-163 



- 38 

- 20 

+ 35 

- 102 

+ 175 

- 65 

- 32 



+ 18 
+ 55 
-137 
+ 277 
— 240 

+ 33 



+ 37 
— 192 

+414 
-517 
+ 273 



From the above table it is seen that the second differences 
are more nearly constant than any of the other series of differ- 
ences. Of equations of form I, 



y = a-\-bx+cx^, 

where x stands for depth and y for velocity, will best represent 
the law connecting the two variables. It should be emphasized, 
however, that the fact that the second differences are nearly 
constant does not show that I is the correct form of equation 
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to be used. It only shows that the equation selected will 
represent fairly well the relation between the two variables. 

It might be suggested that if an equation of form I with 
ten constants were selected these constants could be so deter- 
mined that the ten sets of values given in the table would satisfy 
the equatipn. To determine these constants we would sub- 
stitute in turn each set of values in the selected equation and 
from the ten equations thus formed compute the values of the 
constants. But we would have no assurance that the equation 
so formed would better express the law than the equation of 
the second degree. 

For the purpose of determining the approximate values of 
the constants in the equation 

y = a+bx+cx^ (i) 

from the data given proceed in the following way: 

Let x=X+ico, 

y=¥+yo, 

where xo and yo are any corresponding values of x and y taken 
from the data. The equation becomes 

Y+yo=a+b{X+xo)+c{X+xoy 

=a+bxo+cxo^+(b+2CXo)X+cX^. 

Y = {b+2cxo)X+cX^; (2) 

since yo=a+bxo+cxo^. Dividing (2) by X it becomes 

Y 

— = b+2cxo+cX (3) 

Y 
This represents a straight line when X and — are taken as 

X 

coordinates. The slope of the line is the value of c and the 

intercept the value of b+2cxo. The numerical work is shown 

in the table and the points represented by (x,—\ are seen 
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in Fig. 3. The value of c is found to be —0.76. When a;o=o, 
the intercept, 0.44 is the value of h. For x=X, the value of 



-X 



.4 .5 .6 

Fig. 3. 



yo is taken from the table to be 3.1950, therefore each value of 
Y will be the corresponding value of y diminished by 3.1950. 



X 


y 


X 


Y 


Y 
X 


.441 -.761'' 


a=y-.4^ 


Computed 
y 


.0 


3 1950 





. 0000 




0.0000 


3-1950 


3 - 1948 




I 


3.2299 


I 


0.0349 


0.3490 


0.0364 


3-1935 


3.2312 




2 


3-2532 


2 


0.0582 


0. 2910 


0.0576 


3 


1956 


3-2524 




3 


3.2611 


3 


0.0661 


0. 2203 


0.0636 


3 


1975 


3-2584 




4 


3-2516 


4 


0.0566 


O.1415 


0.0544 


3 


1972 


3.2492 




.■; 


3.2282 


S 


0.0332 


0.0664 


0.0300 


3 


1982 


3-2248 




6 


3.1807 


6 


—0.0143 


— 0.0238 


— . 0096 


3 


1903 


3-1852 




7 


3.1266 


7 


—0.0684 


-0.0977 


— 0.0644 


.3 


1910 


3-1304 




8 


3-0594 


8 


-0.1356 


— 0.1695 


-0.1344 


3 


1938 


3.0604 




9 


2-9759 


9 


— 0.2I9I 


-0.2434 


— 0. 2196 


3 -1955 


2-9752 



10 )31-9476 
0= 3.1948 

The numbers in column 6 were found after b and c were deter- 
mined in Fig. 3. The sum of the numbers in the seventh column 
divided by ten gives the value of a. In the last column are 
written the values of y computed from the formula 



y = T,.ig/\?>-\-.^/\x~']()X^ 



(4) 
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II. y = a+-+-+-+ . . . \ 
Values of — form an arithmetical series and A.''y are constant. 



Another method of determining the constants is illustrated 
in the following example: Let it be required to find an equation 
which shall express approximately the relation between x and y 
having given the corresponding values in the first two columns 
of the table below. 



I 


2 


3 


4 


s 


6 


7 


8 


& 


X 


y 


I 

X 


X 


y 


Ay 


A^y 


2 


Com- 
puted y 


I.O 
1.2 

1-4 
1.6 
1.8 

2.0 
2.2 
2.4 


4.000 
2.889 
2.163 
1.656 
1.284 
1. 000 
0.777 
0.597 


1.0 

0.9 
0.8 
0.7 
0.6 

o-S 
0.4 

0-3 


1. 000 
I. Ill 
1.250 
1.429 
1.667 
2.000 
2.500 
3-333 


4.00 
332 
2.68 
2.07 

i-Si 
1. 00 
0.52 
0.08 


-0.68 
—0.64 

— 0.61 

— 0.56 
-0.51 

— 0.48 
-0.44 


0.04 
0.03 
0.05 
0.05 
0.03 
0.04 


2.00 

i-So 
1. 14 
0.87 
0.67 
0.50 
0.36 
0.25 


4.000 
2.889 
2.163 
1.656 
1.284 
1. 000 
0.777 
OS97 



In colimin 3 are given values of - in arithmetical series 

X 

and the corresponding values of x and y are written in columns 
4 and 5 of the table. The values of y were read from Fig. 4. 
It is seen that the second differences of the values of y given 
in column 7 are nearly constant, and therefore the relation 
between the variables is represented approximately by the 
equation 

This becomes evident if x be replaced by - in I. The law 

X 

may then be stated: 

If two variables, x and y, are so related that when values of - 

CO 

are taken in arithmetical series the n\h differences of the corre- 
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sponding values of y are constant, the law connecting the variables 
is expressed by the equation 



II 



y = a-i 1--5+-5+ 

X XT x^ 






i 














J 




1 


. 


5 


Values of -5 
.6 


7 


. 


3 


. 


) 


1 


V 








































\ 








































3 




\ 




































/ 




\ 


k 






























^ 


/ 




^ 






\ 


s 


























A 









§ 








\ 


S, 






















/ 


















s 


K 
















/ 


-^ 












^ 














^ 


^ 










y 
















1 


















^ 


^ 




k 




































/ 


^ 






































y 


/" 








■~^ 


--~ 


■—. 


^ 
























/ 






















^ 


■ — 



I 



1 1.2 1.4 1,6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 

Values of x 

Fig. 4. 



If in equation (5) - be replaced by X, then 



and 



y = a-^bX-\-cX^, 

' 3;+Aj' = a+&(X+AX)+c(Z+AX)2. 

By subtracting (6) from this equation 

A3/ = tAZ+2c(AX)Z+c(AZ)2; . . 
and from (7) 

Ay + AZy = 6AZ +2c(AX) (Z+AZ) +c(AZ)2. 

Subtracting (7) from (8) 

A2y = 2c(AZ)2; 

L^y 



(6) 



• (7) 
. (8) 



c = 



2(AX)2 
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From column 7 it is seen that the average value of A^.y is 
0.04, and as AX was taken —.1, 



0.04 



Writing the equation in the form 



y-i"H^' 



%'■ 



■2 


y- 


2 
















"71 
























/ 
























/ 






















1 
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/ 
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1 
























1 












.B 
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L 

















.6 ,7 .S .9 1.0 



Fig. 5. 



it is seen that it represents a 
straight line when - andy— 

X 

o 

— are the coordinates. From 

x^ 

Fig. 5 & is found to be 3 
and a to be — i. The for- 
mula is 



,= _,+3,I +,(i^ 



The last column gives the 
values of y computed from 
this equation. 

The following, taken from 
Saxelby's Practical Mathe- 
matics, page 134, gives the 
relation between the poten- 
tial difference V and the cur- 
rent ".4 in the electric arc. 
Length of arc = 2 mm., A is 
given in amperes, V in volts. 



A 

Observed V. 



A 

Computed V . 



1.96 
50-25 



.5102 
50-52 



2.46 
48.70 

■4065 

48.79 



2-97 
47.90 

■3367 
47.62 



3-45 
47-50 

. 2899 

46.84 



3-96 
46.80 

-2525 
46. 22 



4-97 
45-70 

.2012 
45-36 



5-97 
45.00 

■1675 
44.80 



6.97 
44.00 

-1435 
44.40 



7-97 
43.60 

-1255 
44.10 



9.00 
43-50 

.iiii 
43.8s 
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Fig. 6 shows V plotted to — as abscissa. The slope of this 

Ji. 

line is 12.5 divided by .75 or 16.7. The intercept on the V—axis 

is 42. This gives for the relation between V and A 



F=42- 



16.7 



Although the 
points in Fig. 6 do 
not follow the straight 
line very closely the 
agreement between 
the observed and the ^^ 
computed values of V 
is fairly good. 



V 



Fig. 6. 



III. - = a+bx+cx^+do(^+ . . . +oa;". 

y 

Values of x form an arithmetical series and A"- constant. 

y 

If two variables, x and y, are so related that when values of 
X are taken in an arithmetical series the nXh. differences of the cor- 
responding. values of - are constant, the law connecting the variables 

y 

is expressed by the equation 



III 



- = a+bx+cx^+do!^+ . . . -{-ax". 

y 



This becomes evident by replacing y in I by -. The con- 
stants in III may be determined in the same way as they were 
in I. 



IV. y^ = a-\-bx-\-cx^-\-dx^-\- . . . -\-qx". 
Values of X form an arithmetical series and A" y constant. 

If two variables, x and y, are so related that when values of 
X are taken in an arithmetical series the nth. differences of the cor- 
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responding values of y^ are constant, the law connecting the variables 
is expressed by the equation 

rV y^ = a+bx+cx^+dc(?+ . . . +qxr. 

This also becomes evident from I by replacing y by y^. 

The method of obtaining the values of the constants in 
formulas III and IV is similar to that employed in formulas I 
and II and needs no particular discussion. 



CHAPTER II 

V. y = ah\ 

Values of x form an arithmetical series and the values of y a geometrical 
series. 

7/ two variables, x and y, are so related that when values of 
X are taken in an arithmetical series the corresponding values of 
y form a geometrical series, the relation between the variables is 
expressed by the equation 

V y = ab\ 

If the equation be written in the form 

logy=logfl+(log6)a;, 

it is seen at once that if the values of x form an arithmetical 
series the corresponding values of log y will also form an arith- 
metical series, and, hence, the values of y form a geometrical 
series. 

The law expressed by equation V has been called the com- 
pound interest law. If a represents the principal invested, b the 
amount of one dollar for one year, y will represent the amount 
at the end of x years. 

The following example is an illustration under formula V. 

In an experiment to determine the coefficient of friction, /i, 
for a belt passing round a pulley, a load of W lb. was hung 
from one end of the belt, and a pull of P lb. applied to the other 
end in order to raise the weight W. The table below gives cor- 
responding values of a and ix, when a is the angle of contact 
between the belt and pulley measured in radians. 



a 


2 


2ir 
3 


ST 
6 


TT 


6 


4x 
3 


2 


3 


I ITT 

6 


P 


5-62 


6-93 


8.S2 


lo.so 


1 2. go 


15-96 


19.67 


24.24 


29.94 
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The values of a form an arithmetical series and the values 
of P form very nearly a geometrical series, the ratio being 1.23. 
The law connecting the variables is 

P = ab". 

The constants are determined graphically by first writing 
the equation in the form 

log P = log a+a log b 

and plotting the values of a and P on semi-logarithmic paper; 
or, using ordinary cross-section paper and plotting the values 
of a as abscissas and the values of log P as ordinates. Fig. 7 
gives the points so located. The straight line which most 
nearly passes through all of the points has the slope .1733 and 
the intercept .4750. The slope is the value of log b and the 
intercept the value of log a. 



|1.5 

01.0 

m 

V 

^.5 






!4T %'^ %T T %Tr %T %Tr %w ^y^T 
Values o£ a. 

Fig. 7. 



log a =0.4750, 
I'og 6 = 0.1733; 

a = 3, 
6 = 1.49. 



or 



The formula expressing the relation between the variables is 
P = 3(i-49r, 

p = 3g.399a_ 

VI. y = a+b(f. 

Values of x form an arithmetical series and the values of Ay form a 
geometrical series. 

// two variables, x and y, are so related that when values of x 
are taken in an arithmetical series the first differences of the values 
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of y form a geometrical series, the relation between the variables 
is expressed by the equation 

VI y = a-^bc\ 

By the conditions stated the nth value of x will be 

x„ = X]. -\-{n — i) Ax, 

and the series of first differences of the values of y will be 

Ayi, Ayir, Ayxr^, A.y\r^, A.y\r^ . . . Ayir""^. 

The values of y will form the series 

y\, yi+Ayi, y\-\-Ay\-\-rAy],, yi+rAyi+r^Ayi . . . 

yi+Ayi+rAyi+r^Ayi+r^Ayi+ . . . +r"~^Ayi. 

The wth value of y will be represented by 

yn =yi +Ayi . 



From the nth value of x 



« — 1=- 



Ax 



Substituting this value in the above equation there is ob- 
tained 



^yi .^ 



yn=yi+Ayi — 
= a+bc', 

where a stands for yi-\ — —, b for 

i—r i—r 

Let it be required to find the law connecting x and y having 

given the corresponding values in the first two lines of the 

table. 



r ^ , and c for r^^^. 



X 


o 


.1 


. 2 


-3 


-4 


■5 


.6 


-7 


.8 


-9 


I.O 


y 

Ay 

y 


1.300 
0. 140 
1.300 


1.440 

O.IS7 
1-439 


I -597 
0.177 

I -597 


1-774 
0. 200 

1-774 


1-974 
0. 224 

1-973 


2.198 
0-254 
2.198 


2-452 
0.28s 

2-452 


2-737 
0.323 
2.738 


3.060 
0-363 
3 - 059 


3-423 
0.407 

3-421 


3-830 
3-830 
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Since the values of Ay form very nearly a geometrical series 
the relation between the variables is expressed approximately 
by 

y = a+bc'. 

The constants in this formula can be determined graphically 
in either of two ways. First determine a and then subtract 
this value from each of the values of y giving a new relation 

y—a = bc^; 

which may be written in the logarithmic form 

log iy-a) = \og b+x log c, 

and b and c determined as in Fig. 7; or, determine c first and 
plot c* as abscissas to y as ordinate giving the straight line 

y = a+b{c'), 

whose slope is b and whose intercept is a. 

First Method. The determination of a is very simple. 
Select three points P, Q, and R on the curve drawn through 
the points represented by the data such that their abscissas 
form an arithmetical series. Fig. 8 shows the construction. 

P^(xo,a+bc"); 

Q=(xo+Ax, a+bcf'c'^); 

R=(xo+2Ax, a+bc^c"""). 
Select also two more points S and T such that 

S={xo+Ax, a+bc"); 

T=(xo+2Ax, a+b(f'c''^). 
The equation of the line passing through Q and R is 

y= !^ -X !^ -(oco+Ax)+a+bc"c'^. (i) 

Ax Ax 
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y=- 



Ax 



-X — 



Ax 



'(xo+Aa;)+a+6c^°- . (2) 



These Hues intersect in a point whose ordinate is a. For, 
multiplying equation (2) by c'^'^ and subtracting the resulting 
equation from (i) gives 

{i-c^)y = {i-c^)a; 

y = a. 

y 

i. 



R 
'"fa 



.9 1.0 



Fig. 8. 



Fig. 8 gives the value of a equal to 0.2. The formula now 

becomes 

log (y— .2)=log6+«:logc. 

In Fig. 9 log (y — .2) is plotted to .x as abscissa. The slope of 
the line is 0.5185 which is the value of log c, hence c is equal to 
S-S- The intercept is the ordinate of the first point or 0.0414, 
which is the logarithm of b, hence b is equal to i.i. 
The formula is 

3' = o.2+ 1.1(3.3)''. 
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or 



The last line in the table gives the values of y computed from 
this formula. 

Second Method. For any point {x,y) the relation between 
X and y is expressed by 

y=a+bc'', 

and for any other point {x-\-i^x, y+Ay) by 

y-\-^y=a+b(fc'^. 

From these two equations is obtained 

Ay = h(f{c^-i) 

log Ay = log b {c^' — i) +a; log c. 

6 

If now log Ay be 
plotted to X as 
abscissa a straight 
line is obtained 
whose slope is log c. 
The value of c hav- 
ing been determined, 
the relation 



y=a+bi<f) 

will represent a 
straight line pro- 
vided y is plotted 
to c* as abscissa. 




Fig. 9. 



The slope of this line is b and its intercept a. 

VII. log y = a+bc'. 

Values of x form an arithmetical series and the values of A log y form a 
geometrical series. 



// two variables, x and y, are so related that when values of x 
are taken in an arithmetical series the first dijfferences of the cor- 



DETERMINATION OF CONSTANTS 33 

responding values of log y form a geometrical series, the relation 
between the variables is expressed by the equation 

VII log y = a^b(f. 

This is at once evident from VI when y is replaced by log y. 
The only difference in the proof is that instead of the series 
of differences of y the series of differences of log y is taken. 

VIII. y = a^-bx+cd^. 

Values of x form an arithmetical series and the values of A^y form a 
geometrical series. 

If two variables, x and y, are so related that when values of x 
are taken in an arithmetical series the values of the second differ- 
ences of the corresponding values of y form a geometrical series, 
the relation between the variables is expressed by the equation 

VIII y=a+bx+cd\ 
The wth value of x is represented by 

Xn^'Xi + in — ijAx. 
The values of y and the first and second differences may be 



arranged in columns 






yi 


Ayi 




y2 


Ay2 


A2yi 


ya 


Ays 


A2y2 


yi 


Ay4 


A^ya 


ys 


Ays 


A2y4 


ye 


etc. 


etc. 


etc. 







34 EMPIRICAL FORMULAS 

Since the second differences of y are to form a geometrical, 
series they may be written 

The series of first differences will then be 
Ayi, Ayi+A^yi, Ayi+A^yi+rA^yi, Ayi+A^yi+rA^yi+r^A^yi 

Ayi+A^yi+rA''yi+r^A^yi+ . . . +r"~^A'^yi. 

The wth value of y will be equal to the first value plus all 
the first differences. For convenience the rath value of y is 
written in the table below. 

yn=yi 

+Ayi 

+Ayi+A^yi 

+Ayi +A^yi +rA^yi 

+Ayi +A^yi +rA^yi -\-r^A^y\ 

+ Ayi -\-A^yi +rA^yi +r^A^yi +r^A^yi 



+Ayi+A^yi+rA^y\+r^A^yi-{-r^A^yi+ . . . +r'' ^A^yi. 

Adding gives 

-r , i—r^ , i—r^ , i—r^ 



y«=yi + (n — x)Ayi+A^yi\ — 



i—r 



i—r 



I—r J" 



The first two terms on the right-hand side represent the sum 
of all the terms in the first colimin of the value of y„. The 
remaining terms contain the common factor A^yi. The terms 
inside the bracket are easily obtained when it is remembered 
that each line, omitting the first term, in the value of y form a 
geometrical series. It is easily seen that the value of Va may be 
written 
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yn=yi+{n-i)^yi+^^{n-2)-^^{r+r^+r^+ . . . +f^-^) 



i—r i—r 

i—r I—r i—r 



where 



4=3,--m_,5=A3-x+^,andC = 



A^3>i 



(i-r) 



I—r 



A^yi 
(i-rr 



From the value of Xn is obtained 



n—i=- 



Ax 



Substituting this in the value of yn it ia found 

Xg—Xl 

Ax 
= a+bxn+cd^". 

Since x„ and yn stand for any set of corresponding values 
of X and y the resulting formula is 

VIII y = a+bx+c(r. 

In the first two columns of the following table are given 
corresponding values of x and y from which it is required to 
find a formula representing the law connecting them. 



X 


y 


Ay 


A^Ji 


log A^y 


(2.00)* y — I. 


01(2.00)* 


Computed y 


.0 


i.Soo 


.048 


023 


-1-6383 


1. 000 


490 


1.492 


.2 


1-548 


.071 


026 


— I 


5850 


1. 149 


388 


1-550 


■4 


1. 619 


.097 


028 


— I 


5528 


1.320 


286 


1.620 


.6 


1. 716 


.125 


034 


— I 


468s 


I-S17 


184 


I.71S 


.8 


1. 841 


•IS9 


039 


— I 


4089 


1.742 


082 


1. 841 


I.O 


2.000 


.198 


043 


— I 


366S 


2.000 — 


020 


1.999 


1.2 


2.198 


.241 


051 


— I 


2924 


2.300 — 


125 


2.196 


1-4 


2-439 


.292 . 


OS9 


— I 


2291 


2 . 640 — 


227 


2.440 


1.6 


2.731 


•351 


067 


— I 


1739 


3.032 - 


331 


2.735 


1.8 


3-082 
3-500 


.418 








3.482 
4.000 — 


435 


3 -085 


2.0 








S40 


3 506 
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Since the values of x form an arithmetical series and the 
second differences of the values of y form approximately a 
geometrical series, it is evident that the relation between the 
variables is fairly well represented by 

y = a+bx+cd'. 
Taking the second difference 

log A2y = log dd'^-iy+ilog d)x. 



or 



Plotting the logarithms of the second differences of y from 
the table to the values of x, Fig. lo, it is found that log </ = .3000 




.2 Jl ,6 .8 IJD 1.2 lA 1.6 1.8 2S> 
Taliiea o£ <" 

Fig. 10. 

or <i = 1.995, approximately 2. The intercept of this line, 

— 1.6500, is equal to log c{d^ — i)^. 

Since 

d = 2, 

.02239 =£(2"- 1)2, 

C = I.OII. 
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Plotting y- {1.01)2" to x, Fig. lo, the values of a and b 
are found to be 

a= 0.5, 

i= -0.515. 

The formula derived from the data is 

3' = o.5-o.5i5a;+(i.oi)2^ 

In the last colimin of the table the values of y computed 
from the formula are written down. Comparing these values 
with the given values of y it is seen that the formula reproduces 
the values of y to a fair approximation. 

IX. y = io''+'"^+'^'. 
Values of x form an arithmetical series and A^ log y constant. 

If two variables, x and y, are so related that when values of x 
are taken in an arithmetical series the second differences of the 
values of log y are constant, the relation between the variables is 
expressed by the equation 

IX 3, = io«+»^+«=. 

This becomes evident from I when y is replaced by log y. 

log y = a-\-bx-\-cx^, 

which represents a parabola when logy is plotted to x. The 
constants are determined in the same way as they were in 
formula I. 

X. y = ks"g^. 
Values of x form an arithmetical series and values of A' log y form a 
geometrical series. 

If two variables, x and y, are so related that when values of x 
are taken in an arithmetical series the second differences of the 
corresponding values of log y form a geometrical series, the relation 
between the variables is expressed by the equation 

X y=ks'g'\ 
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This becomes evident by taking the logarithms of both sides 
and comparing the equations thus obtained with VIII. X 

becomes 

log y=logk+ (log s)x + (log gjcf. 

This is the same as VIII when y is replaced by log y, a by 
log k, b by log s, and c by log g* 

XL y = 



a+hx+cx^ 
Values of x form an arithmetical series and A"— are constant. 



// two variables, x and y, are' so related that when values oj x 
are taken in an arithmetical series the second differences of the 

corresponding values of - are constant, the relation between the 

y 

variables is expressed by the equation 

X 



XI y=' 



a+bx+cx^ 
Clearing equation XI of fractions and dividing by y 

- = a+bx+coiP. 

y 

This is of the same form as I, and when - is replaced by y 

y 

the law stated above becomes evident. 
If a is zero XI becomes 

I 



b+cx 
which, by clearing of fractions and dividing by y, reduces to 

- = b+cx, 

y 

a special case of III. 

* For an extended discussion of X see Chapter VI of the Institute of 
Actuaries' Text Book by George King. 
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If c is zero XI becomes a special case of XVI, or 

- = a-\-ox, 
J' 

which IS a straight line when - is plotted to x. 

y 

Corresponding values of x and y are given in the table below, 
find a formula which mil express approximately the relation 
between them. 







X 


X 


X 






Y X 




Com- 


X 


y 




A— 


A!_ X 




Y 




-2.5^2 








y 


y 


y 






X y' 




puted y 


o 


o.ooo 
















0.000 


, I 


1-333 


0.075 


100 


050 — 


9 


-2.703 


3 


003 


050 


1.329 




2 


I -143 


O.I7S 


ISO 


050 — 


8 


— 2 . 603 


3 


254 


075 


1. 140 




3 


0.923 


0.32s 


200 


050 — 


7 


-^•453 


3 


504 


100 


0.929 




4 


0.762 


0.525 


250 


050 — 


6 


-2.253 


3 


755 


125 


0.760 




S 


0.645 


0.77S 


300 


051 - 


5 


-2.003 


4 


006 


150 


0.644 




6 


0.558 


I -075 


3S1 


049 - 


4 


-1.703 


4 


257 


175 


0.558 




7 


0.491 


1.426 


400 


047 - 


3 


-1-352 


4 


507 


201 


0.491 




8 


0.438 


1.826 


447 


058 - 


2 


-0.952 


4 


760 


226 


0.438 




9 


0.396 


2.273 


505 


040 — 


I 


-0.503 


5 


030 


248 


0-395 




o 


0.360 


2.778 


545 


054 





0.000 








0.360 




I 


0-33I 


3 323 


599 


056 


I 


0.545 


5 


4SO 


298 


0.331 




2 


0.306 


3.922 


655 


051 


2 


1. 144 


5 


720 


332 


0-305 




3 


0.284 


4-577 


706 


035 


3 


1.799 


5 


997 


352 


0.284 




4 


0.265 


5-283 


741 




4 


2. 505 


6 


262 


383 


0.265 




S 


0.249 


6.024 






5 


3.246 


6 


492 


399 


0.249 



The values of x form an arithmetical series and since the 
second differences of - are nearly constant the values of y will 

y 

be fairly well represented by 

X 
y=- 



or 



a-[-ix-\-cx^' 



- = a+bx+cx^- 

y 



This represents a parabola when - is plotted to x. 
Let X=x — i, 



¥ = --2.778. 

y 
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From these equations are obtained 

x=X+j, 

- = 7+2.778. 

y 

The formula becomes 

F+2.778 = a+6(X+i)+c(Z+i)2 

=o+6+c+(&+2c)Z+cZ2. 
Since the new origin lies on the curve 

a+6+c = 2.778, 




S .2 .3 A .5 .& .1 & 3 hH lA la t3 lA 1.5 X 



Fig. II. 
the equation reduces to 

F=(J + 2C)Z+CZ2, 



or 



- = h + 2C+cX. 



This represents a straight line when — is plotted to X. The 

X 

value obtained for c from Fig. 11 is 2.5. The value of h could 

be obtained from the intercept of this line but the approximation 
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will be better by plotting — i.^x^ to x. In this way is obtained 
the line 

— 2.Kx'^ = a+hx. 

y 

From the lower part of Fig. ii the values of a and b are 
found to be 

ffl = .02S, 
& = .2S25. 

Substituting the values of the constants in XI the formula 
becomes 

X 

.025 + .25253;+ 2. 5x^ 

In the last column of the table the values of y computed 
from this equation are given and are seen to agree very well 
with the given values. 



CHAPTER III 



XII. y = ax?. 

Values of x form a geometrical series and the values of y form a 
geometrical series. 

// two variables, x and y, are so related that when the values of 
X are taken in a geometrical series the corresponding values of y 
also form a geometrical series, the relation between the variables is 
expressed by the equation 

XII y = ax\ 

From the conditions stated equations (o) and (&) are obtained. 

Xn=xir''~^, {a) 

yn = yiR"-\ (b) 

where r is the ratio of any value of x to the preceding one 
and R is the ratio of any value of y to the preceding one. 
Taking the logarithm of each member of (a) 

log x„ =logxi + {n — i) log r, 

loga;„-loga;i 

n—i— ■ — . 

log »• • 

Also by substituting this value of n — i in the value of > in 
equation (b), 

log m— loga!i 

y.=yiR '°^' 



log XI / 1 \ log xa 



=yiR '"6^ \R^°^ 
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where 
and 



log xn 



a=yiR '°er 



io^=R^°^\ 



The following data (Bach, Elastizitat und Festigkeit) refer 
to a hollow cast-iron tube subject to a tensile stress; x represents 
the stress in kilogrammes per square centimeter of cross-section 
and y the elongation in terms of -g^ cm. as unit. 



X.... 

y .... 
log a;, 
logy. 
Comp, 

y- 



9-79 
0-33 
o. 
-o.4»i5 

0.324 



20.02 
0.69s 
1-3014 

-0.1580 

0.714 



40.47 
1.530 
I . 6072 
0.1847 

1.541 



60.92 
2.410 
I . 7847 
0.3820 

2.416 



SI. 37 
3.295 
I . 9104 
0.5178 

3 323 



101.82 
4.18s 
2.0078 
0.6217 

4.252 



204.00 
8.960 
2.3096 
0.9523 

9.132 



408.57 
19.490 
2.6113 
I. 



19.600 



Selecting the values of x which form a geometrical series, 
or nearly so, it is seen that the corresponding values of y form 
approximately a geometrical series, and, therefore, the relation 
between the variables is expressed by the equation 



or 



y = ax , 
log y = log a+S log x. 



If now logy be plotted to log a; the value of h will be the 
slope of the line and the intercept will be the value log a. Fig. 
12 gives h = \.\. In computing the slope it must be remembered 
that the horizontal unit is twice as long as the vertical unit. 
The intercept is —1.5800 or 8.4200 — 10, which is equal to 
log 0.0263. The formula is 

^ = .02635;''^ 

The values of y computed firom this equation are written 
in the last line of the table. They agree quite well with the 
observed values. 



44 



EMPIRICAL rOEMULAS 



XIII. y = a+b log x+c log% 
Values of log x form an arithmetical series and A^y constant. 

If two variables, x and y, are so related that when values of 
log X are taken in an arithmetical series the second differences of 
the corresponding values of y are constant the relation between the 
variables is expressed by the equation 



XIII 



y = a-\-b log x-{-c log%. 



This becomes evident from I by replacing x by log x. The 
law can also be stated as follows: If the values of x form a geo- 



1.4 
1.2 
1.0 

.8 

;•« 
|.2 

|o 

-.2 
-.4 
-.6 

-.8 
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.8 .9 1 1.1 1.2 1.3 \A 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 
Values oi log a: 

Fig. 12. 

metrical series and the second differences of the corresponding 
values of y are constant the relation between the variables is 
expressed by the equation 

y = a-\-b log x-\-c log%. 
If c is zero the formula becomes 

y=a-\-b\Qgx, 
which is V with x and y interchanged. 
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Formula XIII represents a parabola when y is plotted to 
log a;. The constants are determined in the same way as the 
constants in I. 

XIV. y = a-\-lx\ 

Values of x form a geometrical series and values of Ay form a geometrical 
series. 

7/ two variables, x and y, are so related that when the values of 
X are taken in a geometrical series the first differences of the cor- 
responding values of y form a geometrical series, the relation between 
the variables is expressed by the equation 

XIV y = a+hx\ 

As in XII the wth value of x is 

oc„=xir"~'^ (c) 

The series of first differences of y may be written 

Ayi, AyiR, AyiR^, AyiR^ . . . AyiR"'^, 

and the values of y are 

yi, yi+Ayi, yi-\-Ayi-\-AyiR, yi+Ayi+AyiR+AyiR^ . . . 

yx-\-Ay\-\-Ay\R-\-AyxR?-\-AyxR^-\- . . . ■\-AyxR''-^. 

That is the «th value of y will be 

yn = yx^Ayx^AyxR+AyxR^+AyxR^-{- ■ ■ ■ +A3'ii2"-^ 

I— i?""' 
=3/i+Ayi id) 

I — K 

Taking the logarithm of each member of (c), 
log Xn = log iKi + (w — i) log r 
log^^-log^ 

W I — 

log r 
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Substituting this value of w-i in the wth value of y given 
in {£), 



log xn —log Til 



3'n=3'l+A3'l 

Ly 



I-R '°8' 



=yi- 



i-R 

log XI I 1 \loggn 
Ay JJ log r I 2Jlog r J 



I-i? I -22 

= a+6it;„''. 

Let it be required to find the law connecting x and y having 
given the values in the first two lines of the table. 



* 


2 


3 


4 


s 


6 


7 


8 


y 


4.21 


S-2S 


6.40 


7.6s 


8.96 


10.36 


II. 81 


log a; 


.3010 


■4771 


.6021 


.6990 


.7782 


.8451 


.9031 


X 


2 


2.S 


3-I2S 


3.906 


4.883 


6.104 


7.630 


y 


4.210 


4.720 


5.388 


6.290 


7. SIS 


9. no 


11.275 


log a; 


.3010 


• 3979 


.4948 


■S9i8 


.6887 


.7856 




^y 


•Sio 


.668 


.902 


1.225 


I 395 


2.165 




logAji 


-.2924 


-.1752 


— .0448 


.0881 


.2028 


■3358 




y-2.72 


1.49 


2-53 


3-68 


4-93 


6.24 


7.64 


9.09 


log(y-2.72) 


.1732 


.4031- 


.5658 


.6928 


.7952 


.8831 


.9586 


Computed y 


4.21 


S-2S 


6.41 


7.6s 


8.98 


10.36 


II. 81 



In the fourth line values of x are given in a geometrical 
series with the ratio 1.25. In the fifth line are given the cor- 
responding values of y read from Fig. 13. The first differences 
of the values of y are written in the seventh line. These differ- 
ences form very nearly a geometrical series with the ratio 1.336. 
Siace the ratio is nearly constant the law connecting a; and y 
is fairly well represented by the equation 

y = a-\-bx°. 

There are two methods which m^y be employed for deter- 
mining the values of the constants, either one of which may 
serve as a check on the other. 
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The equation of the 
hne passing through P 
and Q is 



First Method. Select three points, A, P, and Q on the 
curve, Fig. 13, such that their abscissas form a geometrical 
series and two other 
points, R and S, such 
that R has the same ^^ 
ordinate as A and the 
same abscissa as P, S k 
the same ordinate as P 
and the same abscissa '■ 
as Q. The points may 
be represented as fol- »' 
lows: » 



A = (xo, a+hxo'); 
P={xor, a+bxoY); 
Q={xor^a+bxoY'); 
R=(pcor, a+bxo'); 
S=(xor^a+bxoV). 
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^ 
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X 










/- 















4 5 

Values of x 







x+a — 



Fig. 13. 

bocoYir'-r) 



xor{r — i) r- 

The equation of the line passing through the points R and 



5 is 



bxo\r'-i) ^ bxo\r'-r) 
y= ) ^x-\-a- 



XQr{r—\) 

These two lines intersect in a point whose ordinate is a. In 
Fig. 13 xo is taken equal to 2 and r equal to 2. The value of 
a is found to be 2.72. The formula then becomes 

y—2.'j2 = bx', 

at 

log (y— 2.72) =log 6+cloga;. 
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In Fig. 14 log (y— 2.72) is plotted to log x and b and c 
determined as in XII. It is seen that the points lie very nearly 
on a straight line. The values of c and b are read from Fig. 14. 

c = i.3; 
log 6 = 9.7840—10; 
b=.6i. 
The law, connecting x and y then is 

y = 2.y2+.6ix^-^. 



.5 

A 
.3 
.2 

bo 

S 

■3 

tu 

3 

t3 

-4 
-:& 



1 \ A. 



A .5 £ .1 .8 .9 1.0 
Values of log x 



1.0 
.9 



I 
.6 a, 



Fig. 14. 

The values of y computed from this formula are written in 
the last line of the table. 

Second Method. From the equation 



y = a+bx' 
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we have 

y+Ay^a+bx'r"; 

Ay = bx'{r''-i) ; 

log Ay = log bir'-i) +c log X. 

This is the equation of a straight line when log Ay is plotted 
to logo;. Fig. 14 shows the points so plotted and from the 
line drawn through them the values of b and c are obtained. 

c = i.3, 
6 = .6i. 

a is found by taking the average of all the values obtained 
from the equation 

a=y—. 6i^-^. 

a is equal to 2.72. 

XV. y = aio''^'. 

Values of x form a geometrical series and A log y form a geometrical 
series. 

If two variables, x and y, are so related that when values of 
X are taken in a geometrical series the first differences of the cor- 
responding values of log y form a geometrical series, the relation 
between the variables is expressed by the equation 

XV y = aio'''\ 

This equation written in the logarithmic form is 

logy=loga+bx''. 

Comparing this with XIV it is evident that if the values of x 
form a geometrical series the first differences of the corre- 
sponding values of log y also form a geometrical series. 

In an experiment to determine the upward pressure of water 
seeping through sand a tank in the form shown in Fig. 15 was 
filled with sand of a given porosity and a constant head of 
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water of four feet maintained.* The water was allowed to flow 
freely from the tank at A. The height of the column of water 

in each glass tube, 
six inches apart, 
was measured. In 
the table below x 
represents the dis- 
tance of the tube 
from the water head 
feet, and y the 



in 

height of the column 
of water in the tube, 
also in feet. It is 
required to find the 
law coimecting x 
Fig. 15. and y. 




Tube 


I 


2 


3 


4 


S 


6 


7 


8 


9 


X 

logy 

X 

logy 

A logy 

log* 

log (-A log y) 

Sx" 

log (.y-bx') 

Computed y 



2.31 
.3636 

.0000 
.3636 
2.314 


5 
'■30 

.3617 

.5 
i.30 

■ 3617 

— .0193 

— .3010 
-1. 7144 

— .0036 
.3653 

2.295 


I 
2 

I 
2 

2 




20 

3424 



20 

3424 

1223 

0000 

9126 

0228 

3652 

19s 


1-5 
2.00 
.3010 

- .0673 
.3683 
1.982 


2.0 
1.66 

.2201 
2.0 
1.66 

.2201 

- .7729 
.3010 

- .1119 

- .1449 
■ 3650 

1.658 


2.5 
1.24 
.0934 

— .2627 

.3561 

A. 264 


3.0 
0.84 

- .0757 

- .4272 
• 3SIS 


3.5 

- .644s 
.3769 
• 525 


0.28 

- .3528 

*°„ 
.28 

- .5528 

- .9201 
.3673 

.278 



In the fifth line values of x are selected in a geometrical 
series and the corresponding values of y written in the next 
line. In Fig. 16 log (—A log y) is plotted to log*. The 
points lie on a straight line. On account of the small number 
of points used in the test we select formula XV on trial. 
From the formula 



it follows that 



y = aio" 
log 3; = logo +6*" 



* Coleman's Thesis, University of Michigan. 
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iog y,c=log a+bxt' 
log yt+i =log a+bxtY 
A\ogyt = bxt''(r''-i) 
log (A log y) = log b(r'-i) +c log x. 

If A log y is negative b is negative, in which case it is only 
necessary to divide the equation by — i before taking .the 
logarithms of the two members of the equation. 




tI .1 

Values oflog x 

Fig. i6. 



The last equation above represents a straight line when 
log (a log y) is plotted to log x. The slope gives the value of 
c and the intercept gives log b(r°-i). From Fig. i6 values 
of b and c are readily obtained. 



C = 2f, 



b= —.02282. 



In the next to the last hne the value of a is computed for 
each value of x from the equation 

log fl = log 3'+ .0228*^'. 
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The average of these values of a gives 

a = 2.314. 
The formula obtained is 

y = (2.3i4)io-''228-''. 

The values of y computed from this equation are written 
in the last line of the table. The agreement is not a bad one. 



CHAPTER IV 
XVI. {x^a){y+b)=c. 

x—xt,, I lie on a straight line. 

y-ykj 
If two variables, x and y, are so related that the points repre- 

(X Xlr\ 
x—Xk, I lie on a straight line, the relation between 
y-yJ 

the variables is expressed by the equation 

XVI {x+a){y^b)=c. 

Let' x—Xk=X, 

y-yt=Y, 

where Xn and y* are any two corresponding values of x and y. 
From the above equations 

X = X+X,r, 

y = ¥+yt. 

Substituting these values of x and y in equation XVI we 
have 

iX+x^+a)(Y+y,+b)=c, 
or 

XY+(y,+b)X+(x^+a)Y+(xH+a)(yt+b)=c. 

Since (xt, yt) is a point on the curve 

{xt+a){yt-\-b)=c, 
and 

XY+{y,+b)X+{x,+a)Y-Q, 

J 

63 
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Dividing the last equation by Y 



X+ (yk+b)Y+Xk+a=o, 



or 



X 

y'' 



I _„_ X +a 



X 



y/c+b" y^+b 
This represents a straight line when Z is plotted to ^. 
The theorem is proved directly as follows: If the points 



x — x , ■ 



'y-yK 



lie on a straight line its equation will be 



X — Xi_ 

y-yk 



■■p{x-x,)+q. 



Clearing of fractions 

x-Xk = p{x-X:){y-yi)-\-q{y-yk)- 

This is plainly of the form 

{x-\-a){y+b)=c. 

The following tables of values is taken from Ex. i8, page 138 
of Saxelby's Practical Mathematics. It represents the results 
of experiments to find the relation between the potential differ- 
ence V and the current A in the electric arc. The length of 
the arc was 3 mm. 



A (am- 




















peres) 


1.96 


2.46 


2-97 


3.45 


3.96 


4-97 


5-97 


6-97 


7. 97 


V (volts 


67.00 


62.7s 


59. 75 


58.50 


56.00 


53.50 


52.00 


51-40 


50.60 


X 





0.50 


Z.OI 


1.49 


2.00 


3.01 


4.01 


5-01 


6.01 


Y 
X 

Y 





-4.2s 


-7.2s 


-8.50 


— 11.00 


-13.50 


— 15.00 


-15.60 


—16.40 




— .1176 


- .1393 


- .1752 


- .1817 


- .2228 


— . 2670 


— .3210 


- .3665 


Com- 




















puted V 


66.99 


62.74 


59-80 


57.80 


56.19 


53-94 


52.44 


51.36 


50-55 



Let A be taken as abscissa and V as ordinate and transfer 
the origin to the point (1.96, 67.00) by the substitution 

X=A-i.g6, 

7=7-67.00. 
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The values of X and F are given in the third and fourth 
lines of the table. The values of — are plotted to X in Fig. 17 





X|>-rl 

ID 
0} 

2-3 

a 
> . 



3 1 

YaluesoJ X 

Fig. 17, 



and are seen to lie nearly on a straight line. It is therefore 
concluded that the formula is 

{V+b)(A+a)=c. 

By the equations of substitution this becomes 

{X+i.g6+a)(V+67.oo+b) =c, 
or 

XF+(67.oo+6)Z+(i.96+6)F=o. 

Dividing by F(67.oo-|-&) 



X 

y" 



X- 



1.96+a 



67.00+6 67.00+& 



The slope of this line is — 

— ^"^ . From Fig. 17 
67.00+0 



67.00+& 



and the intercept is 



Solving these equations 



J^rojn formula 



67.00+& 

1.96+a _ 
67.C0+6 



= ■045; 
095- 



a=o.isi, 
J =-44.78, 

c = 46.89. 
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These values give 

(4+o.isi)(F-44-78).=46.89. 

In the last line of the table are written the values of V com- 
puted from the above formula 

_b_ 

XVIfl. 3/ = aio^+^ 

/ I y y \ 

Points represented by log —, log — lie on a straight line. 

\x-xt yt yt/ 

If two variables, x and y, are so related that the points repre- 
sented by ( log — , log — I lie on a straight line, the relation 

\x-x^ yt yt/ 

between the variables is expressed by the equation 

h 

XVIo y=aid'-^'. 

By the condition stated 

y I 'V 

log — =m log— +J, 

y* x—xt. yt 

where x^ and y* represent any two corresponding values of 
X and y. m is the slope of the line and b its intercept. Clear- 
ing the equation of fractions 

(log y -log yk) {x -Xt) = w(log y — log y*) +b{x-Xt), 
or 

log y{x—Xk-m) = (6+log yt)x-log yt(xt+m) -bxt. 



, ^ (b+log yt)x—log yk(xt+m) —bxt 



x—Xt—m 

^ Ax+B 

x+C 

X+C 

=log »+—;-■ 
x+C 
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Therefore 

b 
-y = Io"'B<'io*+'' 

1) 

= aio"+". 

For the purpose of determining the constants the equation 
is written in the form 

logy=logaH — — , 
x+c 

(logy-log a)ix+c) =b, 
Let 

logy=logF+logyt, 
and 

x=X+Xk. 
Then follows 

(log F+log yt-log a){X+Xt+c) = b, 

Xlog F+log F(a;t+c)+X(log3't-loga)+(log)'i-log a)(xi+c) =b. 

But 

(log yic — log a) (act+c) = b, 

since the point {xu, yk) lies on the curve. 

X log F+log F(xt+c) +Z(log yi-log a) =o. 
Dividing this equation by X 

log F= - {xu-\-c) -^- +log a-log yu. 

Replacing log F and X by their values 

log Z = _ (a;t+c)^^ log -+log a-log y^. 

From this it is seen that if log — be plotted to log ^ 

yt x-xt yt 

a straight line is obtained whose slope is —{xt+c) and whose 
intercept is log a— log y^. If the slope of the line is represented 
by M and the intercept by B 

c= —M—Xt, 
loga=B+logyt. 
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By writing XVIa in the logarithmic form 

logy = 6-— — hloga 
x+c 

a line is obtained whose slope is b. 

XVII. y = ae"+be'^. 

(yi+i yt+aX 
> 1 
y* yt I 

lie on a straight line whose slope, M, is positive and intercept, B, is negative, 
and M^+4B positive. 

// iwo variables,- x and y, are so related that when values of »^ 
are taken in an arithmetical series the points represented by 

/yk+iyk+2\ ^^g ^^ ^ straight line whose slope, M, is positive 

\ y^' y" / 

and whose intercept, B, is negative and also M^-\-^B is positive 

the relation between the variables is expressed by the equation 

XVII y = ae'^+be^. 

Let {xk, yi), {x+hx, yk+\), (a;i+2Aa;, yk+2) be three sets of 
corresponding values of x and y where the values of x are taken 
in an arithmetical series. We can then write the three equations, 
provided these values satisfy XVII. 

y, = ae''''^+be'^^\ (i) 

y*+i = ae''=V^^+&e^V^^ (2) 

y,+2 = ae''"'e^'^='+be^=''e^^^'' (3) 

Multiplying (i) by e'^^^ and subtracting the resulting 
equation from (2) 

y,+x-e'^''y, = be^'''{e'^^''-e'''^) (4) 

Multiplying (2) by e'^'^^ and subtracting from (3) 

y^+2-e'^^y,+^=be'^K^^''''-e^c+d)^x). . . (s) 

Multiplying (4) by e'^^^ and subtracting from (5) there 
results 



or 
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y*+2-(e'=^^+e'^^^)y.+i+e(«+*^^>'.=o,. . . (6) 
y«:+2 ^ rgC&x_i_gdAx\ yt^+1 _ g(c +d)&x 

y* yk 

The values of c and d are fixed for any tabulated function 
which can be represented by XVII, and therefore, the last 

equation represents a straight line when ^^^^^ is plotted to ^^^. 

yk yk 

The slope of the line is 

and the intercept is 

5 = _ g(c +d)Ax^ 

It is seen that M is positive, B negative, and M^+4B posi- 
tive, for 

ilf2 = ^CAX_^ ^^{c +d)Ax_^^dAx^ 



and 



45= -4e(«+*^*, 



lf2+45 = (eCAx_grfAa;)2_ 

In the first two lines of the table are given corresponding 
values of x and y. It is desired to find a formula which will 
cfpress tjhe relations between them. 



* 


I.O 


IS 


2.0 


2.5 


3.0 


3-5 


4.0 


4-5 


S.o 


y 


+ .3762 


+.0906 


-.1826 


- .4463 


- .7039 


- .9582 


— i.2iig 


-1.4677 


— 1.7280 


Vk+l 

yk 






+ ■241 


-2.0I5 


+2.444 


+I.S77 


+1.361 


+1.265 


+ 1.2II 


yk+2 
yk 






-.48s 


-4.926 


+3.855 


+2.147 


+1.722 


+I.S32 


+ 1.426 


g-A12x 


+ .662 


+ •539 


+ .439 


+ .359 


+ .290 


+ .236 


+ .192 


+ .IS7 


+ .127 


Computed y 


+ .319 
+ .371 


+ .071 
+ .087 


-.131 
-.185 


- .295 

- .447 


- .429 

- .704 


- :S38 

- .957 


- .626 
—1. 210 


- .698 
-1.464 


- -757 
-1.723 



Plotting the points represented by (^^^, ^^j, Fig. 18, 
a straight line is obtained whose equation is 



yk+2 „„)'*+ 1 „/; 

— = 1-97'^:; -96, 



y^ 



yk 
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Jf = i.97, 
B=- .96. 

Since M is positive, B negative, and M^+^B positive, it 
follows that the relation between the variables is expressed 
approximately by XVII. It has been shown that the slope 
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.3 .4 
Values of e"' 
Fig. 19. 



of the line is equal to e'^^^+e'^^, and the intercept is equal to 

eJ'«+'')= .96. 

From these equations are obtained the values of c and d, 

c=-.247, 
d = .i6s. 
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The formula is now 

Dividing both sides of this equation by c"'* gives the 
equation 

which represents a straight line when ye~-^^^" is plotted to 
g-.4i2i -pjjg values of these quantities taken from the table 
are plotted in Fig. 19 and are seen to lie very nearly on a straight 
line. 

This line has the slope 2.cx) and intercept — i.oi. Sub- 
stituting these values of a and b in the formula it becomes 

It is seen that the errors in the values of y computed from 
this formula are in the third decimal place. The values are as 
good as could be expected from a formula in which the con- 
stants are determined graphically. For a better determination 
of the constants the method of Chapter VI must be employed. 

XVIII. y = e'^ic cos hx-\-d sin hx). 
Values of x form an arithmetical series, and the points ( , ) 

\ y* yt / 

lie on a straight hne. Also 'M.^-\-i,B is negative. 

7/ two variables, x and y, are so related that when values of 
X are taken in an arithmetical series the points represented by 

iHl^ yj^\ lie on a straight line whose slope M and intercept 

yt yic / 

B have such values that M^+^B is negative, the relation between 
the variables is expressed by the equation 

XVIII y=ef"^{c cos bx+d sin bx). 

Let X and yt be any two corresponding values of the variables. 
We have the three equations 

yk=e'"'{c cos bx+d sin bx), (i) 

yt+i =e''V^^[c cos (bx+bAx)+d sin (bx+bAx)] 
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= e''V^''[c(cos bx COS &Ax — sin hx sin SAx) 

+rf(sin hx cos SAx+cos &x sin 6Aa;)] 

= e"'V^^[(c cos iAx+^ sin 6Ax)cos bx 

+ (^ cos SAx — c sin 6Ax)sin 6a;j. . (2) 

The value yK+2 can be written directly from the value of 
yii+ 1 by replacing Ax by 2Ax. 

3't+2 = e''V'"^^[(c cos 2&Ax+</ sin 2&Ax)cos bx 

+ (<i cos 2JAa;— c sin 2&Aa;)sin &x] (3) 

Subtracting (i) multiplied by e'"^^(c cos bAx+d sin &Ax) 
from (2) multiplied by c we have 

cyk+i—e'"^{c cos b^x+d sin.SAx)ys 

_^g<MgaAi^^ cos &Ax — c sin J Ax) sin &x 

—de'-^^'^'^ic cos SAx+i sin &Ax)sin 6x 
= -(c2+rf2)e=v^^sin&Axsinix (4) 

Similarly 

c3'*+2 — e^^^^Cc cos 26AX+J sin 26Ax)3'i: 

Multiplying equation (4) by e"'^'' sin 26AX and subtracting 
it from (s) multiplied by sin 6Ax 

c sin bi^xyic+2 — e'^'^^{c cos 2&Ax sin &Ax+J sin 2&Ax sin 6Ax)ys 

—ce'^" sin 26Ax3't+i +e^'"^''(c cos &Ax sin 2&Ax 

-\-d sin 6Ax sin 26Ax)3;i; = o. 
Simplifying 

c sin 6Ax3't+2— ce""^^ sin 26Axyft+i+ce^'"^ sin JAxyt=o. 

Dividing by c sin SAx^t, 

^^ = 2 cos bAxe'"^'^^-e'"'^\ 
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The values of a and b will be fixed for any tabulated function 
which can be represented by XVIII, and therefore, the last 

equation represents a straight line when ^^ is plotted to 
^^. The slope of the line is 



and the intercept 



M = 2e°^^cos6Aa;, 



It is evident that M^+4B is negative. 

It is possible that in a special case M^+^B might be zero, 
but then b would be zero and hence 



which is formula V. 

Corresponding values of x and y are given in the first two 
columns of the table below. It is required to find a formula 
which will represent approximately the relation between them. 





y 


yk+1 

yk 


yk+i 
yk 


^.08. 


COS bx 


tan bx 


cos 63c 


y 


Com- 


X 


e-o^^cos bx 


puted y 





+ .300 






I . 0000 


i.OOOO 


.0000 


i.OOOO 


+ .300 


+ .308 


I 


+ .Oil 






1.0833 


+ .8646 


+ .5812 


+ .9366 


+ .012 


+ .018 


2 


- .332 


+ .04 


— i.ii 


1.1735 


+ .4950 


+ 1.7556 


+ .5809 


- .571 


- .327 


3 


- .636 


—30.2 


-S7.8 


I. 2712 


— .0087 


-114-59 


— .0111 


+57.3 


- .634 


4 


- .803 


+1.92 


+ 2.42 


I. 3771 


— .5100 


— 1.6864 


- .7023 


+ 1.143 


— .804 


S 


- .761 


+1.26 


+ 1.20 


I. 4918 


- .8732 


■ 5581 


— 1.3026 


+ .584 


- .761 


6 


- .48s 


+ .95 


+ .60 


I.6I6I 


- .9998 


+ .0175 


-1.6159 


+ .300 


- .48s 


7 


- .017 


+ .64 


+ .02 


1.7507 


- .8557 


+ .604S 


— 1. 4981 


+ .oil 


— .012 


8 


+ -537 


+ .04 


— I. II 


1.8965 


- .4797 


+ 1.8291 


— .9098 


- .590 


+ .54s 


9 


+1.027 


-31.6 


—60.4 


2.0544 


+ .0262 


— 38.1880 


+ .0538 


+ 19. 08 


+I.03S 


ID 


+1.298 


+I.91 


+ 2.42 


2.2255 


+ .5250 


— i.6212 


+I.I684 


+ 1 .III 


+ 1.299 



In Fig. 20 the points represented by / Z?±i^ Z^:^ \ are plotted. 

\ y^ yk / 

They lie very nearly on the straight line whose equation is 



^ = 1.875^-1.175. 
jk yk 
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Since (1.878)^—4(1.18) is negative the relation between the 

variables is expressed 
approximately by the 
equation 



y = e'^'^ic cos hx-\-d sin ix). 

It was shown that 
the slope of the line is 
equal to 2 (cos lAx)e''^ 
and the intercept equal 

equal to unity we have 

26" cos 5 = 1.875, 



'^"" = i-i7S, 



log cos b = 9.9370—10, 

6 = 30° 10' ap- 
proximately, 

a = .08. 

The formula is now 
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Fig. 20. 



-y = e<'s^(c cos 3oia--f(f sin 30^0;), 

where 30^ is expressed in degrees. 

Dividing the equation by e""^ cos 30^3; 



-pg; j-5 = c-f-J tan 30^35°, 

e""^ cos 3o|x 



which is a straight line when 



-^—5 is plotted to tan 



e"°^ cos 3oia;' 

3o-6-a;°. In Fig. 21 these points are plotted and are seen to 
lie nearly on a straight line whose slope is — .496 and intercept 
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.308. Two of the points are omitted in the figure on account 
of the magnitude of the coordinates. Substituting the values 
of constants just found in the formula the equation expressing 
the relation between x and y is 

<y==g.08i^_^Qg COS 3oia;°-.496 sin 30^*°). 
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Fig. 21. 



yt+i y*+2 



The last column 
in the table gives 
the values of y 
computed from the 
equation. The 
agreement with the 
original values is 
fairly good. 

In case c is zero 
XVIII becomes the 
equation for damped vibrations, y = de"^ sin bx. 

XIX. y = axl+bx''. 
Values of x form a geometrical series, and the points 

\ Vk Vk 

lie on a straight line, whose slope, M, is positive, and whose intercept 
B, is negative, and M'+4.B positive. 



If two variables, x and y, are so related that when values of 
x are taken in a geometrical series the points represented by 

lyk+iyic+2\ ^^^ ^ ^ straight line whose slope, M, is positive, 
\ yk yt / 

and intercept, B, negative, and also M^+4B positive, the relation 
between the variables is expressed by the equation 

XIX y = ax''+bx''. 

Let X and yt be any two corresponding values of the variables. 
The following equations are evident: 

yt = ax'+bx*, (i) 

yt+i = ax'r'+bx^r', (2) 



66 EMPIRICAL FORMULAS 

3'.-r2 = axV^+6it;V, (3) 

Jk+i-r^y^ ^bx'ir'-r'), (4) 

yk+2-r'yt+i =boifr^ir^-r'') (s) 

Multiplying equation (4) by y" and subtracting it from 
equation (5) there results 

yt+2—r'yt+i—r^y!c+i+r''+''yt=o, 
or 

y*+2 ^ (r---\.r^) y^fc+i yc+g^ 
y« yt 

It is seen that the slope of this line is positive and the inter- 
cept negative, and M^+^B positive. 

In the table* below, the values of x and y from a; = .05 to 
a; = .5 5 are taken from Peddle's Construction of Graphical 
Charts. 











yic+i 


y!c+i 


.55 


85 


y 


Com- 


X 


y 


X 


y 


yt 


yk 


X 


X 


x'' 


puted y 


■ OS 


.283 


OS 


283, 






.192 


078 


1.470 


-283 




10 


.402 


10 


402 






.282 


141 


1.426 


.402 




IS 


.488 










■352 


199 


1-385 


.488 




20 


.556 


20 


556. 


1.420 


' 1-965 


-413 


255 


1-347 


-556 




25 


.613 










.466 


308 


I-315 


.612 




30 


.6^8 










-SI6 


359 


1.276 


.658 




3S 


•69s 










-561 


410 


1-238 


.697 




40 


•730 


40 


730 


'i'383' 


1. 816 


.609 


4S9 


1.208 


•730 




45 


• 7S7 










-64s 


507 


1. 174 


• 757 




SO 


.780 










.683 


SSS 


1. 142 


.780 




55 


.800 










.720 


602 


1 .114 


•799 




60 


.814 










•755 


648 


1.078 


.814 




65 


.826 










-789 


693 


1.047 


.826 




70 


•83s 










.822 


738 


1. 016 


.835 




75 


.840 










.8S4 


783 


0.984 


.840 




80 


• 84s 


80 


84s 


1-313 


1 . 520 


-88s 


829 


0-9S5 


.846 



In column 3 the values of x are selected in geometrical ratio 
and the corresponding values of y are given in column 4. The 

points 1^^^^^^^) are plotted in Fig. 22, and although the 

\ yt yt / 

* See Rateau's " Flow of Steam Through Nozzles." 
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three points do not lie exactly on a straight line the approxima- 
tion is good. The slope of the line is 4.10 and the intercept 
which ffive the eniiations 



—3.86 which give the equations 



2'^+ 2" =4.10, 



-.c+d _ 



=3-86. 
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2.00 



1.90 



1.70 jS 



1.6a 



1.50 



1.30 1.32 1.34 1.36 1.38 1.40 1.12 

Values of ^ 

Fig. 22 AND Fig. 23. 

Solving these equations the values of c and d are found to be 
6 = 1.40, 

'^ = -55- 



The formula now is 



y=ax'-*+6r^^ 
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Dividing both members of this equation by a;'^ 

which represents a straight line when -^ is plotted to t^. 

The slope of the line is equal to a and the intercept equal to b. 

From Fig. 23 

a=- .685, 

6 = 1.522. 

The formula after the constants have been replaced by their 
numerical values is 

y = i.$22x-^^ — .68$x^'^''. 

The last column of the table shows that the fit is quite 
good. 

If the errors of observation are so small that the values 
of the dependent variable can be relied on to the last figure 
derivatives may be made use af to advantage in evaluating the 
constants in empirical formulas. But when the values can not 
be so relied on, or when the data must first be leveled graphically 
or otherwise, the employment of derivatives may lead to very 
erroneous results. This will be illustrated by two examples 
worked out in detail. 

The first step in the process is to write the differential equa- 
tion of the formula used and then from this equation find the 
values of the constants. 

Consider the formula 

y = e°^(c cos bx+d sin bx). 

Looking upon a and b as known constants and c and d 
as constants of integration, the corresponding differential 

equation is 

y" — 2ay' + (a^+b^)y = o. 

Dividing this equation by y 

^ = 2a5^-(a2+62) 

y y 
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which, if the data can be represented by XVIII, represents a 
straight line whose slope is 2a and whose intercept is —{a^+V^). 
Corresponding values of x and y are given in the table. 





y 


y' 


y" 


t. 


?! 




t 


gMx 


cos .08a: 


tan -08a: 


y 


X 






6-081 COS .o8j 










y 


y 


De- 


Min- 




















grees 


utes 













-■ .3000 

















1. 0000 


1. 0000 


.0000 


-f -3000 


X 


+ .2750 










4 


35. 02 


I. 0618 


.9968 


.0802 


+ -2598 


1 


+ .2441 


-.0342 


-!cio68 


— .1401 


— .0279 


9 


10.04 


I. 1275 


.9872 


.1614 


+ -2193 


3 


+ .2065 


-.0439 


-.0066 


- .1976 


-•0319 


13 


45. o5 


I. 1972 


.9713 


.2447 


+ .1776 


4 


-f .1622 


— . 0481 


— .0078 


- .296s 


-.0481 


18 


20.08 


I. 2712 


.9492 


.3314 


+ -1344 


S 


+ .1102 


-.05S7 


-.0075 


- .5054 


-.0681 


22 


55- 10 


1-3499 


.9211 


.4228 


-f -0886 


6 


+ .0506 


-•063s 


-.0086 


-1.2549 


— .1700 


27 


30.12 


1-4333 


.8870 


.5206 


+ -0398 


7 


- .0175 


— .0721 


-.0080 


+4.1200 


--•4571 


32 


05.14 


1 .5220 


.8472 


.6270 


— -OI36 


8 


- .0937 


-.0805 


-.0087 


-- .8591 


--.0928 


36 


40.16 


I.6161 


.8021 


.7446 


- .0723 


9 


- .1786 


— .0894 


— .0091 


- - .5011 


--.0510 


41 


IS. 18 


I. 7160 


-7519 


.8771 


- -1384 


10 


— .2726 


— .0985 


— .0091 


-- .3247 


--.0334 


45 


50.20 


I. 8221 


.6967 


1.0296 


- -2147 


II 


~ -^w 


-.1078 


-.0085 


-- .2869 


--.0226 


50 


25.22 


1.9348 


-6372 


I . 2097 


- -3047 


12 


— .4881 


-.1168 


— .0087 


-- .2393 


--.0178 


55 


00.24 


2.0S44 


-5735 


1.4284 


- .4143 


13 


- .6093 


-.1257 


-.0091 


- - . 2063 


--.0149 


59 


35.26 


2.181S 


-5062 


1.7036 


- -5518 


14 


- .7396 


-•1348 


— .0089 


-- .1823 


--.0120 


64 


10.28 


2.3164 


-4357 


2.0659 


- .7328 


IS 


- .8788 


-.1435 


-.0084 


+ .1633 


--.0096 


68 


45.30 


2.4596 


-3624 


2.5722 


- -9859 


16 


—1.0264 










73 


20.32 


2.6117 


.2867 


3.3414 


-1.3707 


17 


-1.1814 










77 


55-34 


2.7732 


.2098 


4-6735 


—2-0306 



The values of y' and y" are obtained by the formulas 

y'n = — 7 (yn-2-8y„-i+ 8y„+i-y„+2), 
i2n 

y"n = 71^ {yn-2 — i6y„-i +soy„ — l6yn+l+yn+2), 

I2fr 

where h = Ax — i. These formulas are derived in Chapter VI. 

Plotting the points represented by (—.—], Fig. 24, it is 

seen that they lie nearly on a straight line whose slope is .12 
and intercept —.01. Therefore 



We have then 



2a = .i2, 

a^+P = .oi, 

a = .06, 

& = .o8. 

y=e-°^''{c COS .o8x+d sin .08a;). 
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Dividing this equation by e"^* cos .08a; 

=c+d tan .08*. 



e cos .o8x 
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Fig. 24. 



This represents a straight line when 



e"^^ cos .08a; 



is plotted 



to tan .08a;. The slope is d and intercept c. From Fig. 25 

c= .3. 
The law connecting the variables is represented by 
3/ = 6°°'' (.3 cos .08a;— .5 sin .o8x). 
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The values of y computed from this formula agree with 
those given in the second column of the table. 
Consider formula XIX 

The corresponding differential equation is 





t/ 
























e0.6ar 
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Va 


I 


I tan 


3 

.08a; 




■1 









s 



.4 
2 
O" 
-.2 o 






:6 
-iS 
■1..9 
■1.2 
-1.4 
-1.6 
-1.8 
-3.0 



Fig. 25. 

where c and d are known constants and a and h constants of 
integration. The differential equation represents a straight line 

9, ff t 

when ^^ is plotted to — . The slope is c-\-d--L and the inter- 
X y 

cept is —cd. 

The values of x and y in the table below are the same as 
those given in the discussion of XIX. 
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xy £ 


y^ 


X 


y 


y' 


y" 


y 


y 


■OS 


283 












lO 


402 












IS 


488 


1-503 


-6.933 


.462 — 


320 




20 


556 


1.240 


-4-133 


.446 - 


297 




25 


613 


±.015 


-4.967 


.414 - 


506 




30 


658 


0.803 


-3-267 


-366 


447 




35 


69s 


0.720 


—0.400 


•363 - 


071 




40 


730 


0.623 


-3-533 


-341 - 


774 




45 


757 


0.492 


-1.500 


.292 


401 




50 


780 


0.433 


— 1.067 


-277 - 


342 




55 


800 


0.338 


-2.633 


.232 — 


996 




60 


814 


0-2SS 


-0.633 


.188 - 


280 




6S 


826 


0.213 


— 1.200 


.167 


614 




70 


835 


0135 


-1.767 


.113 -I 


037 




75 


840 












80 


84s 











The values of y' and y" were computed by means of the 
formulas used in the preceding table. In Fig. 26 the points 

—, ) are plotted and, as is seen, the points do not deter- 

y y I 

mine a line. It is clear that the constants can not be determined 
by this method. 

XlXfl. y=dy^if. 
Points represented by (*„, log 5!^M ije on a straight line. 

7/ two variables, x and y, are so related that when values of x 
are taken in a geometrical series the corresponding values of y are 

such that the points represented by (x„, log^^ j lie on a straight 
line, the relation between the variables is expressed by the equation 
XlXa 3; = aa;V. 

Using logarithms: 

log ^B =log a-\-b log Xn-\-x„ log c, 
log yn+i =log a+b log x„+rxn log c+b log r. 
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Subtracting the first equation from the second 
log 5l=±i = {r-i)xn log c-\-h log r. 

By plotting log 2!2±i to a:„ a line is obtained whose slope is 
yn 

(r — i)log c, and since r is known c can be determined. 











o 




-.1 

-.2 
-.3 
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-.4 
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= -.5 
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g-6 






















t7 






























o 




-;8 






















-.9 






















-1.0 


























° 









Values of 
Fig. 26. 



XV ' 



From the first equation 

log yn-x„ log c = b log a;„+log a. 

If then log jn—Xa log c be plotted to log x„ a line is obtained 
whose slope is b and whose intercept is log a. 



CHAPTER V 

XX. y=0o+0i COS «+a2Cos 2a;+a3 003 3*+ . . . -^-arCosrx 

+bisinx+bism 2x+b3smsx+ . . . +6r_isin (y— i)*. 
Values of y periodic. 

The right-hand member of XX is called a Fourier Series 
when the number of terms is infinite. In the apphcation of the 
formula the practical problem is to obtain a Fourier Series, 
of a hmited number of terms, which will represent to a sufl&ciently 
close approximation a given set of data. The values of y are 
given as the ordinates on a curve or the ordinates of isolated 
points. 

In what follows it is assumed that the values of y are periodic 
and that the period is known. 

We will determine the constants in the equation , 

y = ao+ai cos x+a2 cos 2X+a3 cos ^x,+bi sin x+b2 sin 2X, 

so that the curve represented by it passes through the points 
given by the values in the table. 



X 


o° 


60° 


120° 


180° 


240° 


300° 


360° 


y 


I.O 


1-7 


2.0 


1.8 


i-S 


0.9 


1.0 



Substituting these values in the equation we have the fol- 
lowing six linear relations from which the values of the six 
constants can be determined : 

i.o = ao+ 0,1+ 02+03, 

1-7 =ao+|ai- -202-03-1 -bi-] ^62, 

2 2 

2.0 = 00-|01-|02 + 03-| ^61 ^&2, 

2 2 
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1.8 = ao— ai+ a2—az, 

i.5 = ao-^ai-|o2+a3 -h-\ — -1)2, 



2 2 



^A._^3, 



0.9 = 00+501 — 202—03 =^61 -b2- 

2 2 

Multiplying each of the above equations by the coefficient 
of oo (in this case unity) in that equation and adding the result- 
ing equations we obtain (i) below. Multiplying each equation 
by the coefficient of oi in that equation and adding we obtain 
(2). Proceeding in this manner with each of the constants 
a new set of six equations is obtained. 

600 = 8.9 (i) 

301 = -1.25 (2) 

302= -.25 (3) 

603 = . 10 (4) 

3&i = .6sv'3. . (5) 

3^2 = . 15^3 (6) 

O0=ft, fll=— A-, 02=— T^, 03=CT, 

The equation sought is 

3' = !^- A cos x- 1*2- cos 2x+-h cos 3a;+M'^3 sin a;+^V3 sin 2X. 

It reproduces exactly each one of the six given values. 

The solution of a large number of equations becomes tedious 
and the probabihty of error is great. It is, therefore, very 
desirable to have a short and convenient method for com- 
puting the nimaerical values of the coefficients.* 

*The scheme here used is based upon the 12-ordinate scheme of 
Runge. For a fuller discussion see "A Course in Fourier's Analysis and 
Periodogram Analysis " by Carse and Shearer. 
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Take the table of six sets of values 



X 


o° 


60° 


120° 


180° 


240° 


300° 


y 


yo 


yi 


y2 


ys 


y4 


y5 



where the period is 2ir. 

For the determination of the coefficients the following six 
equations are obtained: 

3'o = ao+ ai+ ^2+03, 

yi=ao+iai-ia2-a3-\ — -bi-\ — ^62, 

2 2 

^2 = ao - iai ^ 1^2 +33 -\ — -bi -b2, 

2 2 

yz=ao — ai + ^2 — 03, _ 

yi = ao-^ai-^a2+a3 -bi-\ — ^62, 

2 2 

y5 = ao+iai—^a2 — a3 ^61 ^62. 

2 2 

Proceeding in the same way as was done with numerical 
equations the following relations are obtained: 

6ao=yo+ yi+ y2+y3+ 3/4+ 3)5, 

301=3/0+ bi- iy2-y3- lyi+hys, 
3(12= yo- iyi- hz+ys- hi-hys, 
6a3=yo- yi+ >'2-3'3+ yi— ys, 

^~3 



3bi 



3^2 



-\ — -yi-{ — -y2 -^yi- 

2 2 2 



-ys, 



, '^3 ^3 

+^yi ^y2 

2 2 



2 2 



(a) 



For convenience in computation the values of y are arranged 
according to the following scheme: 



yo 
ys 



yi 
yi 



y2 



Sum 
Difference 



Wo 



Vi 
Wi 



■02 
■W2 
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Vl 

V2 



Wo 



Wl 

W2 



Sum pa 

Difference 






Sum 
Difference 



ro 



Sl 



6ao=po+pi, 
3ai=ro+isi, 
2,a2=po-lpi, 
6az = ro-Si, 



3f> 



_v^ 



■ri, 



(i) 



3^2 = -^qi. 



It is evident that the equations in set (b) are the same as 
those in set (a). 

For the numerical example the arrangement would be as 
follows: 

i.o 1.7 2.0 

1.8 1.5 0.9 







Vo 


2.8 




3-2 


2.9 










Wo 


-.8 




.2 


I.I 








2.8 




3-2 






— 


8 


.2 




2.8 




2.9 
6.1 




ro 






I.I 


Po 


— 


.8 


1-3 


qi 






•3 
6ao = 

3^1 = 

3?2 = 

6a3 = 
361 = 
3^2 = 


= +8 
= — I 

= + 

= + 
= + 


Sl 

.90, 

■25, 

■25, 
.10, 

•65^/3, 

•IS^/•3■ 






-•9 
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It is seen that the computation is made comparatively simple. 
The values of the v's are indicated by vo, the first one. The 
values of the ^'s, etc., are indicated in the same way. 

8-OEDiNATE SCHEME. The formula for eight ordinates which 
lends itself to easy computation is 

y = ao+ai cos d+a2 cos 26+03 cos 30+04 cos 46 
+61 sin 5+62 sin ad+bs sin ^8. 

For determining the values of the constants eight equations 
are written from the table : 






45° 


90° 


135^^ 


180° 


225° 


270° 


31S" 


yo 


yi 


y2 


ys 


y* 


ys 


ye 


3'7 



yo = ao+ ai+(i2+ as+ai, 

V 2 V 2 '\^2 V2 

yi = ao-\ ai as-a^-l bi+b2-\ is, 

2222 

3/2 = 00 — 02 +ai+ bi — bs, 

V2 , V2 , V2, , , V2 , 

y3 = ao fli H a3 — a4:-\ oi — 02H 203, 

2 2 2 2 

y4 = ao— 01+02— 03+04, 

"^2 , V2 V2, , . \/2, 

3'5 = Oo oi H 03 — 04 01+62 63, 

2 2 .2 2 

3'6 = Oo -02 +04-5l + 63, 

,, , V2 V2 V2, , ^2, 

y7=ao-\ oi 03 — 04 Oi — 02 03. 

2222 

From which are obtained the following eight equations: 

8oo=yo+ yi+y2+ y3+yi+ ys+y6+ yi, 

, V2 V2 V2 V2 
4ai=yo-\ yi >'3-3'4 ys H yv, 

2 2 '2 2 

402= yo -3/2 +3'4 -ye,' 

'^2 , -\/2 , V'2 V2 
4«3=yo — — yi i ys—yi-i ys y?. 
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804 =^0- yi+y2- ya+yi- ys+ye- yi, 



461= — yi+y2-\ ys 

2 2 

462 = 

463 = 



V2 V2 

— --ys-ye y?, 

2 2 



yi - ys + ys - y?, 



V2 , V2 

— yi—y2-i y3 

2 2 



V2 , V2 

y5 +y6 y7. 

2 2 



For the purpose of computation the values of y will be 
arranged as follows: 

yo y\ y2 yz 

yi 3'5 y& yi 



Sum 
Difference 



Siun 
Difference 



I'D 



I'D 
1>2 



111 
Wi 

1>3 



V2 
W2 



W3 



Wo 



po pi Sum 

qo qi Difference 

8ao = po+ pi, 

V-2 

4ai=ro-{ si, 

2 

4a2 = qo, 

V~2 

4a3=ro 51, 



ro 



8ai=po—pi, 

V2 

4h = r2-\ ri, 

2 

4h = qi, 

V~2 

463=— '•2 H n. 

2 

The process will be made clear by an example 

B o 45° 90° 135° 180° 225° 270 



Wl 
W3 



ri 

Sl 



° 315° 



W2 



r2 



360° 



4 —2 —I 



— I 
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For computation the arrangement is as follows: 
4—2—12 
3 3-12 





■vo 


7 


I 


— 2 


4 






Wo 


I 


-5 












' 7 


I 




I 


-5 







— 2 
5 


4 
5 


ro 









^ 


I 


-5 





?o 


9 


-3 


Sl 

&ao = io, 

4^1 = 1— fv^2, 

4^2 = 9, 

4a3 = i+fV2, 

8a4 = 0, 

4&1 = -IV2, 

4^2= -3, 

4J3=-fV2. 




-s 





The formula becomes 

)' = i.25 — .634 cos 0+2.25 cos 20+1. 134 cos 30 
— .884 sin 9— .75 sin 20— .884 sin 30. 

io-Ordinate Scheme 





yo 


yi 


y2 


3'3 


yi 










ys 


ys 


yr 


ye 


3'5 




Sum 


Vo 


Vi 


V2 


VI3 


Vi 


vs 




Difference 




Wl 


W2 


W3 


Wi 








■Vo 


Vl 


V2 






Wl 


W2 




V5 


n 


V3 


Sum 




Wi 


Wz 


Sum 


Po 


Pi 


P2 


h 


h 


Difference 


go 


?i 


?2 


Difference 


w:i 


OT2 
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ioao=po+pi+p2, 
5(ii = qo+Ciqi+C2q2, 
5a2=po+C2pi - Cip2, 
5<is=qo-C2qi-Ciq2, 
5a4:=po — Cipi+C2p2, 

ioa5 = qo-qi+q2, 
5bi=Sili+S2l2, 
Sb2=S2mi+Sim2, 
Sh=S2h—Sil2, 
$bi= SiMi— S2m2- 



In the above equations 

Ci = cos36°, 
C2 = cos 72°, 



5i = sin36°, 
' 52 = sin 72°. 



In the schemes that follow, as in the lo-ordinate one, only 
the results will be given. 



12-Ordinate Scheme 







yo 


yi 


y2 


3'3 


3'4 




y5 


Je 








yii 


yio 


y^ 


y» 




yi 




Sum 


Vo 


Vl 


V2 


Vz 


Vi 




va 


% 


Difference 






Wl 


■W2 


Wz 


Wi 




Ws 






vo 


Vl V2 


V3 




Wl 


W2 


W3 








V6 


V5 Vi 






•W5 


Wi 









Sum po pi p2 p3 

Difference qo qi 92 



Simi ri r2 

Difference si S2 



rs 



po 

p2 



p3 



ri 



qo 
92 



Sum 



lo 



Difference 



t2 
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J2ao = lo + h, 

6ai=qo-\ — -qi+M2i 

2 

6a2=po-p3+i(pi-p2), 

6fl3=fc, 

6ai=po+p3-i{pi+p2), 
6o5=?o — -qi+iq2) 

2 

I2a6=^~^i) 
6bi = ^ri-\ — ^r2+r3, 

2 

6b2=—^{si+S2), 

' 2 

663=^1, 

6h= — ^(51-52), 
2 

2 



16-ORDINATE Scheme 

Jo yi ^2 ^3 yi 75 ys yv ys 

yi5 yii y:3 712 yn ym ys 

Sum I'D Vi V2 Vz Vi V5 1'6 V7 V8 

Difference wi W2 wz Wi W5 wq wr 









V2 
V6 


VZ 

V5 


Vi 


Sum 
Difference 


pa 
90 


pi 
11 


P2 
?2 


pz 

93 


P4^ 
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Wl W2 Wz Wi 

W7 Wq W5 



Sum ri r2 rs n 

Difference ^i 52 sz 

po pi p2 lo h 

pi p3 h 



Sum h h h Sum to 

Difference mo mi Difference xq 

x(>ao=to-\-h, 

■\/~2 

8ffli =g'oH q2+Ciqi+C2q3, 

V~2 

8a2 = woH mi, 

2 

8(13 = Qo g'2 — Ciqz -\-C2q1, 

2 

8ai=xo, 

8a5 = ?o ?2 +Ciq3— C2qi , 

8a6=mo mi, 

2 

8a7 = qo-\ q2 — Ciqi—C2q3, 

2 

i6a8=io—h, 

V2 

861 = r4:-\ r2 +C\r3 +C2ri, 

2 

862 = 52H (51+53), 

2 

V2 

8&3= —r4:-\ r2+Ciri—C2r3, 

854 = 51—53, 
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V~2 

?>h = ri r2+Ciri —C2r3, 

2 

8&6= -52 H (^i+^s), 

2 
V2 

8b7=—ri r2 + Cirs + C2ri . 

2 

Ci = cos 22|° = sin 67!°, 
C2 = sin 22|° = cos 67!°. 



20-Oedinate Scheme 

yo yi y2 ys yi ys ye yi y?, y^ yic 
yw yi8 yn yie yu yu y\z yi2 yn 



Slim 


vo 


Vl 


V2 V3 


Vi 


V5 


V6 


V7 Vs 


Vg Z)io 


Difference 




Wl 


W2 W3 


\ Wi 


W5 


Wo 


Wl Ws Wg 






vo 


Vl 


V2 


vz 


Vi 


v& 








Vio 


va 


Vs 


Vl 


Vo 






Sum 


po 


Pi 


P2 


pz 


pi 


P5 




Difference 




go 


gi 


12 


qz 


?4 










Wl 


W2 


Wz 




Wi 


W5 








■Wg 


Ws, 


Wl 




W6 






Sum 




ri 


r2 


rs 




ri 


Th 




Difference 




Si 


S2 


sz 




Si 








po 


pi 


p2 






?o 


gi 


?2 




P5 


pi 


Ps 






qi 


?3 




Sum 


h 


h 


h 


Sum 


ko 


hi 


k2 


Difference 


mo 


mi 


W2 














lo 






Wo 


mi 




ri n 




h 






W2 








rs 




h 
















Sum 


to 




Sum 


no 


m 


Sum 


Ol 03 
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Si 


S2 


S4 


S3 


gl 


g2 


h 


fe 



Sum 

Difference 
2000 = to, 

ioai = qo+qi sin 72°+g2 sin S4°+?3 sin 36°+54 sin i8°, 
ioa2=mo+mi sin 54°+W2 sin i8°, 
T.oa3=qo—q3 sin 72° — 94 sin 54°+5'i sin 36°— 52 sin 18°, 
ioa4 = fo— ^2 sin 54°+/i sin 18°, 
1005 = ^0 — fe, 

ioa6=wo— OT2 sin 54°— mi sin 18°, 
ioa7 = qo+q3 sin 72° — 94 sin 54° — 51 sin 36° — 52 sin 18°, 
iofl8 = fc— ^1 sin S4°+/2 sin 18°, 

ioaQ=qo—qi sin 72°+g2 sin 54°— 93 sin 36°+g4 sin 18°, 
2oaio=no—ni, 

lobi^r^+ri sin 72°+/3 sin S4°+r2 sin36°+>'i sin 18°, 
1062 =g2 sin 72°+gi sin 36°, 

1063 = — r5+>'2 sin 72°+?'i sin 54°— ^4 sin 36°+y3 sin 18°, 
iobi = hi sin 72°+fe sin 36°, 
1065 = 01—03, 

io&6=gi sin 72°— g2 sin 36°, 

1067= —r5—r2 sin 72°+^ sin SA°+i'i sin 36°+r3 sin 18°, 
1063= —>^2 sin ']2°+hi sin 36°, 
i6bQ = r5—ri sin 72°+r3 sin 54°— >'2 sin 36°+ri sin 18°. 

24-OEDiNATE Scheme 

yo yi y2 yz yi ys y& yi y?, y^ yw yn yi2 
y23 3'22 ^21 y2o yi9 yi8 yn ym yis yu yis 

Sum vo vi V2 V3 Vi V5 Vq V7 vs vo vm Vn 1)12 

Difference %vi ws ws W4 W5 we wj wa wg wio wn 



86 



EMPIRICAL FORMULAS 



Vo Vi 
Dl2 Vn 



V2 

no 



V3 

V9 



Vi 



V7 



Sum 
Difference 



Sum 
Difference 



Sum 
Difference 



po 



Wl 



pi p2 

qi ?2 

ze'2 

Wio 



p3 

qz 






pi 
qi 

Wi 
■W8 



P5 

W7 



p6 



We 



ri 

Sl 



r2 

S2 



Si 



ri 

Si 



S5 



po pi p2 p3 

p6 pS pi 

lo h h h 

mo mi m2 



re 

Sl S2 
S5 Si 



lo 
h 



h 
k 



Sum 
Difference 



go 
ho 



gi 
hi 



Sum 
Difference 



Sum 
Difference 



Wo 

W2 



kl k2 
til »2 

mi 



eo 
fo 



ei 



24ao=go+gi, 
i2ai=qo+iqi+iV2q3+^V;^q2+Ciqi+C2q5, 

1202 = Wo + |ot2+i V3OT1, 

i2a3=qo-qi+iV2{qi-q3-q5), 
i2ai = ho+\hi, 

I2a5 = ?0 + C2?l-jV3g2-|V2y3+i?4 + Cl?5, 
1206 =/o, 

I2a7 = go-C2?l-|V3^2 + |Vi?3 + |?4-Cig5, 

i2a8 = ?o-igi, 

I2a9=go-g4 + iV2(-yi+j3+g5), 
12^10 = Wo + |W2-|V3Wl, 

I2ail=?0-Cl?l+iV3g2-|V2g3 + ^?4-C2g5, 

24ai2 = ki-hi} 



Si 



kz 
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i2&i=C2ri+|r2+|v'2r3+|V3r4+Cir5+/'6, 
i2b3 = r2-r6+iV2{ri+ri—r5), 

I2&4 = |V3(mi+«2), 

i2b5 = Ciri+^r2-iV2r3-iV;^n+C2r5+r6, 

12hQ = ki — kz, 

i2b7 = Ciri-^r2-iV2r3+W3n+C2r5-r6, 

I268=fV3(wi-W2), 

i2bQ=r6-r2+i^2(ri+r3-r5), 

I2bio = S3+l{si+S5)—^^3(s2+S4), 

i2bn = C2ri—ir2+^V2r3—h^3U+Cir5-rQ, 

2V2 

C2= ^~^ = .25882. 
2V2 

As an illustration let it be required to find a Fourier series 
of 24 terms to fit the data given in the table below. 



x" 


y 


x" 


y 


X' 


y 


x" 


y 


00 


149 


90 


159 


180 


178 


270 


179 


IS 


137 


105 


178 


19s 


170 


28s 


i8s 


30 


128 


120 


189 


210 


177 


300 


182 


4S 


126 


135 


191 


225 


183 


31s 


176 


60 


128 


ISO 


189 


240 


181 


330 


166 


• 75 


135 


i6S 


187 


255 


179 


345- 


160 



149 137 128 126 128 135 159 178 189 191 189 187 178 

160 166 176 182 185 179 179 181 183 177 170 

vo 149297 294 302 310 320 338 357370374366357178 

?2'i -23 -38 -50 -54 -50 -20 -I 8 8 12 17 
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149 297 294 302 310 320 338 

178 357 366 374 370 357 

po 327 654 660 676 680 677 338 

§'0 —29 —60 —72 —72 —60 —37 

-23 -38 -50 -54 -50 -20 

17 12 8 8-1 

n — 6 —26 —42 —46 —51 —20 

si -40 -50 -58 -62 -49 

327 654 660 676 -40 - 50 -58 

338 677 680 -49 - 62 



h 665 1331 1340 676 ki —89 —112 -58 

mo —II —23 —20 n\ 9 12 

665 1331 -II -23 

676 1340 —20 



go 1341 2671 eo -31 -23 

h -II -9 /o 9 

The formula becomes 

3; = 167. 167 — 19.983 cos a;— 3. 410 cos 2a;+s.47o cos 3a; 

— 1.292 cos 40;+. 249 cos sa;+.75 cos 6a;-f-.3io cos 7a; 
+.458 cos 8a;— .304 cos 9a; — .090 cos loa; — .243 cos iia; 

— .083 cos 12a;— 12.779 sin a; — 16.624 sin 2a; — .323 sin 3a; 
+ 1.516 sin 4a;+i.46i sin 531-2.583 sin6a;+.32i sin 7a; 

— .216 sin 8a;+.676 sin 9a;— .459 sin loa;— .639 sin iia;. 

In what precedes the period was taken as 2-a. This is not 
necessary; it may be any multiple of 27r. The process of finding 
a Fourier series of a limited number of terms which represent 
data whose period is not 27r will be best set forth by an example. 
In the table below the period is x/3 and the values of y are 
given at intervals of -rfi?,. The 12-ordinate scheme can be 
used by first making the substitution 

x=\d or = 3a;. 
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x" 


0° 


y 


x° 


0° 


y 


a:° 


e° 


y 


oo 

lO 

20 

30 


00 
30 
60 
90 


+27.2 

+34-5 
+ 21. S 
+ 10. 1 


40 

5° 
60 
70 


120 
ISO 
180 
210 


+9.8 

+8.S 
+0.2 

-7.1 


, 80 

90 

100 

110 


240 
270 
300 
330 


-ii-S 
-I7-S 

-17.2 

+ i.S 





27.2 


34-5 


21-5 


10 


I 9 


.8 8.5 0.2 






1-5 - 


-17.2 - 


17 


5 -II 


•5 -7-1 


I'D 


27.2 


36.0 


4-3 - 


-7.4 -I 


.7 1.4 0.2 


Wl 




33-0 


387 


27 


6 21 


•3 iS-6 






27.2 


36.0 




4-3 


-7-4 




po 


0.2 


1.4 




-1-7 






27.4 


37-4 




2.6 


-7-4 




go 


27.0 


34-6 




6.0 






33-0 


38-7 


27.6 






27-4 37-4 




15-6 


21.3 








2.6 -7.4 


ri 


48.6 


60.0 


27.6 




/o 


30.0 30.0 


5l 


17-4 


17-4 








48.6 27.0 
27.6 6.0 



21.0 



21.0 



The formula is 



;y = 5+9.994 cos 0+8.7 cos 2^+3.5 cos 39+.006 cos 50 
+ 17.31 sin e+s.023 sin 20+3.5 sin 30-. 01 sin 50. 

Replacing 6 by its value 3a;, 

y = S+9-994 cos 3X+8.7 cos 6X+3.5 cos 9X+.006 cos iSic 
+17.31 sin 3a;+s.o23 sin 6^+3.5 sin 9a;-. 01 sin 15a;. 



CHAPTER VI 

EMPIRICAL FORMULAS DEDUCED BY THE METHOD 
OF LEAST SQUARES 

In the preceding chapters we computed approximately the 
values of constants in empirical formulas. The methods em- 
ployed were almost wholly graphical, and although the results 
so obtained are satisfactory for most observational data, other 
methods must be employed when dealing with data of greater 
precision. 

It is not the purpose of this chapter to develop the method 
of least squares, but only to show how to apply the method) to 
observation equations so as to obtain the best values of the 
constants. For a discussion of the subject recourse must be 
had to one of the numerous books dealing with the method of 
least squares.* 

It was found in Chapter I that the equation 

y=a+bx+cx^ (i) 

represents to a close approximation the relation between the 
values of x and y given by the data 



X 


y 


X 


y 


o 


3-1950 


■5 


3.2282 


.1 


3.2299 


.6 


3.1807 


.2 


3-2532 


•7 


3.1266 


■3 


3.26II 


.8 


3-0594 


4 


3-2516 


•9 


2-9759 



* Three well-known books are: Merriman, Method of Least Squares; 
Johnson, Theory of Errors and Method of Least Squares; Comstock, 
Method of Least Squares. 
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where x represents distance below the surface and y represents 
velocity in feet per second. 

Substituting the above values of x and y in (i), the following 
ten linear observation equations are found : 

c+ 0&+ oc = 3.i95o, 
a+.ib+.oic =3.2299, 
a+.2&+.o4c=3.2S32, 
a+-3^+-09C = 3-26ii, 
a+.4J+.i6c = 3.2516, 
a+.56+.25c=3.2282, 
a+.6&+.36c=3.i8o7, 
ff+-7^+-49C = 3-i266, 
, a +.86+. 64c =3.0594, 
a+.9J+.8ic = 2.9759. 

Here is presented the problem of the solution of a set of 
simultaneous equations in which the number of equations is 
greater than the niunber of unknown quantities. Any set of 
three equations selected from the ten will suffice for finding 
values of the unknown quantities. But the values so found 
will not satisfy any of the remaining seven equations. Since all 
of the equations are entitled to an equal amount of confidence 
it would manifestly be wrong to disregard or throw out any one 
of the equations. Any solution of the above set must include 
each one of the equations. 

The problem is to combine the ten equations so as to obtain 
three equations which wiU yield the most probable values of 
the three unknown quantities a, b, and c. It is shown in works 
on the method of least squares that the first of such a set of 
equations is obtained by multiplying each one of the ten equa- 
tions by the coefl&cient of a in that equation and adding the result- 
ing equations. The second is obtained by multiplying each one 
of the ten equations by the coefiicient of J in that equation 
and adding the equations so obtained. The third is obtained 
by miiltiplying each of the ten equations by the coefficient of 
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c in that equation and adding the equation so obtained. The 
process of computing the coefficients in the three equations is 
shown in the table. The coefScients of a, b, and c are represented 
hy A, B, and C respectively, and the right-hand members are 
designated by N. The number S, which stands for the numeri- 
cal sum of ^, J5, C and N, is introduced as a check on the work. 
It must be remembered that this method of finding the values 
of the constants holds only for linear equations. 

The sum of the numbers in the colimin headed A A = 2.4^4 
= io. The sum of the numbers in the column headed AB = 
2^5=4.5. The sum of the numbers in the colxmin headed 
AC = 'SAC = 2.85. Also the sum of the numbers in the colimm 
headed ^iV = S^iV = 3i. 7616. These sums give the coefficients 
in the first equation.* The second and third equations are 
obtained in the same way. 

The three equations from which we obtain the most probable 
values of the constants are : 

10 o -I- 4.56 -f-2.85c =31.7616; 
4.5a -1-2.856 +2.02SC =14.08957; 
2.8sa -f2.0256-fi.S333c= 8.828813. 

These are called normal equations. From them are obtained 

«= +3-19513; 
b = + .44254; 
c=- •76531- 
I'he check for the first equation is 

i:AA + i:AB+-ZAC+i:AN = I,AS=4g.jxi6; 
for the second equation 

1:AB + -2BB + -2BC+2BN = I:BS = 25.46457; 
for the third equation 

S^C-f 25C-F2CC -I- 2CiV = 2C5= 15.237113. 
* Cf. Wright "and Hayford, Adjustment of Observations. 
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AA 



AB 



AC 



o 

.01 
.04 
.09 
.16 

•2S 
•36 

■49 
.64 
.81 



AN 



1950 
2299 

2532 
2611 
2516 
2282 
1807 
1266 
OS94 



2 ■9759. 



AS 



4-1950 
4-3399 
4-4932 
4.6511 
4.8116 
4.9782 

5 - 1407 
5-3166 

5-4994 
5-6859 



4-5 



2.8s 



31.7616 



49.1116 



AB 


BB 


BC 


BN 


BS 


















.01 


.001 


.32299 


■43399 




.04 


.008 


.65064 


.89864 




.09 


.027 


■97833 


I ■ 39533 




.16 


.064 


1.30064 


1.92464 




•25 


-I2S 


1.61410 


2.48910 




•36 


.216 


I . 90842 


3.08442 




-49 


•343 


2.18862 


3.72162 




.64 


■512 


2.44752 


4^ 39952 




.81 


.729 


2.67831 


S^ii73i 


4-5 


2-85 


2.025 


14.08957 


23 ■ 46457 


AC 


BC 


CC 


CN 


CS 




' 















.0001 


.032299 


■043399 






.0016 


.130128 


.179728 






.0081 


• 293499 


.418599 






■ 0256 


.520256 


.769856 






■ 0625 


. 807050 


I . 244550 






.1296 


I. 14505 2 


1.850652 






■ . 2401 


1.532034 


2.605134 






.4096 


I. 958016 


3.519616 






.6561 


2.410479 


4^ 605579 


2.8s 


2 -02s 


I 5333 


8.828813 


15^ 237113 
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The formula is 

y =3. 19513+442543;-. 76531^:^. 

For the purpose of comparison the observed values and the 
computed values are written in the table, v (called residual) 
stands for the observed value minus the value computed from 
the formula. 





Observed 


Ccmpute 1 






X 


y 


y 


V 


»2 





3-1950 


3.1951 


— .0001 


. OOOOOOOI 


.1 


3 • 2299 


3.2317 


-.0018 


.00000324 




2 


3-2532 


3-2530 


+ , 0002 


. 00000004 




3 


3.261I 


3-2590 


+ .0021 


.00000441 




4 


32516 


3-2497 


+ .0019 


.00000361 




S 


3.2282 


3-2251 


+ .0031 


.00000961 




6 


3.1807 


3-1851 


-.0044 


.00001936 




7 


3.1266 


3.1299 


--0033 


. 00001089 




8 


3-0594 


3-0594 


.0000 


. 00000000 




9 


2.9759 


2-9735 


+ .0024 


.00000576 



+ oooi .00005493 

This method derives its name from the fact that the sum of 
the squares of the residuals is a minimum. A discussion of 
this will be found in the books referred to above. 

In case the formula selected to express the relation between 
the variables is not linear the method of least squares cannot 
be applied directly. In order to apply the method the formula 
must be expanded by means of Taylor's Theorem. Even when 
the formula is linear in the constants it may be advantageous 
to make use of Taylor's Theorem. In order to make this trans- 
formation clear we will apply it to the formula just considered. 

Suppose that there have been found approximate values of 
a, b, and c, ao, &o and co, say, then it is e\'ident that corrections 
must be added in order to obtain the most probable values of 
the constants. Let the corrections be represented by Aa, Ab, 
and Ac. And let 

a = ao+Aa, 

b-=bo+Ab, 

c = co+Ac. 
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The formula was 

y=a+bx+a(^. 
This may be written 

y=fia, b, c)=f(ao+Aa, bo+Ab, co+Ac). 

Expanding the right-hand member 

fiao+Aa, bo+Ab, co+Ac)=f{ao,bo,co)+^Aa+^Ab+^Ai 

dao doo dco 

+pF|J^(Aa?+|^,(A6)2+|^^(A.)2+-9|-(AM6) 
[2 ld(kr db(f dc(f doodbo 

+^^{AaAc)+-f^—iAbAc)]+ . . . 
9ao9co doodco J 

where -^ stands for the value of the partial derivative of 
doo 

{{a, b, c) with respect to a and ao substituted for a, — ^ stands 

for the value of the second partial derivative of J{a, b, c) with 
respect to a and ao substituted for a, etc. If ao, bo, and co have 
been found to a sufficiently close approximation the second and 
higher powers of the corrections may be neglected. 

^-; 

3/ 2 
9co 

The formula becomes 

y-fiao, bo, co)^-^Aa+^[Ab+^Ac, 
d(k> oOo oco 



or 

y—icio+ box + cox^) =Aa+ xAb + x^Ac. 
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Selecting for the values of oo, bo, and co those found in 
Chapter I, the new set of observation equations are 

Aa+ oAS+ oAc= .0002, 

Afl+.iA6 + .oiAc= — .0013, 

Affl+.2A& + .04Ac= .0008, 

Ao+.3A&+.09Ac= .0027, 

Aa+.4A& + .i6Ac= .0024, 

A(i+.SA6 + .25Ac= .0034, 

Aa+.6A&+.36Ac= -.0045, 

Afl+.7A6+-49^c= -.0038, 

Ao+.8AS + .64Ac= —.0010, 

Ao + .9A6+.8iAc= .0007. 
From these are obtained the three normal equations 
ioAa+4.5 A&+2.85 Ac=-.ooo4, 
4.5AaH-2.85 A&+2.025 Ac=— .00203, 

2.8sAa+2.02SA6+i.5333Ac= -.002059. 

Solving 

-Aa = +.00033, 

A6 = +.00254, 
Ac =-.00531, 
which added to the values of ao, bo, and co, give 

«= 3-i9Si3> 
b= .44254, 

^=-•76531- 
the same as just found. 

The above process may be applied to linear equations con- 
taining more than three constants. But as the method of pro- 
cedure is quite evident from the above the general statement 
of the process will be made with reference to equations con- 
taining only three constants. 
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Let the observation equations be represented by 
aix+biy+ciz=ni pi, 
a2X+b2y+C2Z=n2 p2, 
a3X+b3y+C3Z=n3 p3, 



amXr\-bmy + CmZ=nm pm- 

The normal equations will then be 

S^a^ • X + "Lpab ■y+Xpac-z = "Zpan, 
Xpab ■ X + Xpb^ ■ y + llpbc ■ z = 'Zpbn, 
Xpac ■x+'2pbc-y+ lipc^ ■ z = 'Zpcn, 

where a, b, c, and n are observed quantities, and x, y, and z 
are to be determined, pi, p2, ps ■ • ■ pm are the weights assigned 
to the observation equations. In the problem treated at the 
beginning of the chapter the weight of each equation was taken 
as unity. 

It was stated on a preceding page that when a formula to be 
fitted to a set of observations is not linear in the constants it 
must be expanded by Taylor's Theorem. 

Take as an illustration a problem considered in Chapter IV. 

The formula considered was 



y= 


=f(A, B, 


m, n) 


^Ax'^+Bx", 




3/. 
dAo 


=it;'"°, 






9/. 


=x'", 






3/. 
dmo 


=Aox"" 


' log X, 




df . 
dno 


=Box'" 


log a;; 
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y =f{A ,B,m,n) =f(Ao, Bo, m, no) +^A^ +^^^ 

-L ^■^-Aot+— Aw; 



dmo dno 

y-Mo, Bo, mo, m) =^AA+^AB+^Am+^An. 
dAo dBo onto o^o 

The observation equations will be of the form 

dAo dBo dmo dno 

Assume the approximate values found in Chapter IV. 
A= 1.522, 
5= -.68s, 

m= -SS, 
n= 1.4. 



,.55 

ll4 



X 

Ao 

Bo 

logo; 

Agx^'logx.. 
Box"" log I. 



.OS 


.10 


.15 


.20 


.25 


.19 


.28 


• 3S 


.41 


•47 


.02 


.04 


.07 


.10 


.14 


1.522 










- .68s 










— 2.996 


-2.303 


-1.897 


— 1 . 609 


-1.386 


- .88 


- -99 


— 1.02 


— 1. 01 


- 98 


■03 


.06 


.09 


.12 


■14 



■30 
•52 
.19 



— 1 . 204 

- .94 

.15 



..14 



X 

log a; 

Aox"^logx.. 
Box'" log X. 



■35 


.40 


•45 


•56 


.60 


.64 


•23 


.28 


•33 


-1.050 


— 0.916 


-0.799 


- .90 


- .84 


- .78 


.16 


•17 


.18 



■5° 
.68 
.38 
-0.693 
- .72 
.18 



55 
72 

43 

598 

66 



The new observation equations become 

.19A.4+.02A5— .88Aot+.o3Aw= .0004, 
.28A.4-I-.04A5— .99A?»+.o6Aw= .0002, 
.3SA^+.o7A5— i.02Am-l-.09Aw= —.0001, 
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.4iA^ + .ioA5 — i.oiAw+.i2A«= .0000, 
.47A^+.i4A5— .98Am+.i4Aw= .0013, 
.S2A^+.i9A5— .94Aw+.isA«= — .0001, 
.S6AA+.2^AB— .9oAw+.i6Aw= — .0019, 
.6oA-4+.28A5— .84Am+.i7A«= — .0016, 
.64A4+.33A5— .78Aw+.i8Aw= — .0001, 
.68A^ + .38A5— .72Aot+.i8Aw= — .0001. 
.72A^+.43A5— .66Am + .iS>An= .0011. 

From these the four normal equations are obtained 

2.96oA^ + i.32iAB — 4.637AOT+ .8o6A« = —.00071, 

1.321A.4+ .642A5 — 1.802AOT+ .359Aw= —.00031, 

— 4.637A^ — i.8o2A5+8.737Am— 1.253AW = +.00085, 

.806A.4+ .359A5— 1.253AW+ .22iA«= —.00023. 

From which 

AA = —.0068, 

A5 = +.oii2, 

Am= —.0022, 

An= —.0070. 

These corrections being applied the final formula becomes 
y = i.Si52r5478_.6738a;i-393. 



CHAPTER VII 

INTERPOLATION.— DIFFERENTIATION OF TABULATED 
FUNCTIONS 

Interpolation 
In Chapter II we found that the formula 

XL y= * 



02S+.2S2sa:+2.sa;2 



represents to a fair degree of approximation the values of y 
given by the data. Any other value of y, within the range of 
values given, can be obtained in the same way. This rests on 
the assimiption that the formula derived expresses the law con- 
necting X and y. For example, the value of y corresponding 
to a; = i.o5 will be 

I. OS 
\ = =0.^4.';. 

•" .025 + .2S2S(l.OS) + 2.S(l.05)2 ^^^ 

When a formula is used for the purpose of obtaining values 
of y, within the range of the data given it is called an inter- 
polation formula. Interpolation denotes the process of calcu- 
lating under some assumed law, any term of a series from values 
of any other terms supposed given.* It is evident that empirical 
formulas cannot safely be used for obtaining values outside 
of the range of the data from which they were derived. 

* For a more extended discussion of the subject the reader is referred 
to Text-book of the Institute of Actuaries, part II (ist ed. 1887, 2nd ed. 
1902), p. 434; Encyklopadie der Mathematischen Wissenchaften, Vol. I, 
pp. 799-820; Encyclopedia Britannica; T. N. Thiele, Interpolationsrechnung. 

As to relative accuracy of different formulas, see Proceedings London 
Mathematical Society (2) Vol. IV., p. 320. 
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There are two convenient formulas for interpolation which 
will be developed.* 

The first one of these requires the expression for yx+n in 
terms of y^ and its successive differences, y^ represents the 
value of a function of x for any chosen value of x, and yx+n 
represents the value of that function when x+n has been sub- 
stituted for x. 

yx+i=yx+Ayx; 

yx+2 =yx+Ayx+A(yx+Ayx) 

=yx+2Ayx+A^yx; 
yx+3=yx+2Ayx+A^yx+A(yx+2Ayx+A^yx) 

=yx+3Ayx+3A^yx+A^yx; 
yx+i=yx+3Ayx+3A^yx+A^yx+A(yx+3Ayx+3A^yx+A^yx) 

= yx+4Ayx + 6A^yx +4A^yx + A^yx. 

These results suggest, by their resemblance to the binomial 
expression, the general formula 

yx+n=yx-\-nAyx-\ — ^n --A^yx^ — ^ f^ ^ASy.+etc. 

1^ \3 

If we suppose this theorem true for a particular value of n, 
then for the next greater value we have 

, ^ ,w(w — 1)^„ , «(w — i)(w— 2)^, , , 

yx+n+\=yx■\■nAyx^ — ^-1 '-A^yx^ — ^ ^ ^A^^/^+etc, 

P \3 

+Ayx+nA^yx+-^. — ^A^y^+etc, 

=yx+{n+i)Ayx+ ^ A^yx+- V^ ^A^y^+etc 

P |3 

The form of the last result shows that the theorem remains 
true for the next greater value of n, and therefore for the next 

* See Chapter III, Boole's Finite Differences. 
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greater value. But it is true when n=4, therefore it is true 
when w = 5. Since it is true for » = 5 it is true when n = 6, etc. 
If now o is substituted for x and x for n, it follows that 



yx=yo+xAyo- 



^(^-i)a2„ I x(x-i)(x-2) 



-A^yo 



\3 



A^yo+etc. 



If A"y.=o, the right-hand member of the above equation 
is a rational integral function of x of degree n — i. The formula 
becomes 

, . ,x(x—l)^, , x{x — l){x — 2) ^^ , 

y^ = yo+xAyo+^ — -A^yo+ ■ ,, -A^yo + 



x{x—l)(x — 2) 






(i) 



Formula (i) will now be applied to problems. It must not 
be forgotten that in applying this formula x is taken to represent 
the distance of the term required from the first term in the series, 
the common distance of the terms given being taken as unity. 

I. Required to find the value of y corresponding to a; = .42 5 
having given the values under XIX. In the interpolation 
formula ^ = .5. 

yo y\ y2 yz 



Aya. . 
A^yo. 
^yo- 

y=yfi-\- l^yo - { A^yo + r^ a^^-o 

= .730+.oi35 + .ooo5-|-.oooi 

= .744. 

This is the same as given by XIX. 

2. Find the value of y corresponding to a; = 2. 3 
formula will have the value | if we take yo=—. 

X = 2. 



■ 730 


•757 


.780 


.800 


.027 


.023 


.020 




— .004 


-.003 






.001 









X in the 
1826 when 
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yo 


yi 


3*2 


yz 


yi 


3'5 


-.1826 


- 4463 


-•7039 


-.9582 


— 1.2119 


-1.4677 


-.2637 


-.2576 


-•2543 


-•2537 


- -2558 




.0061 


■0033 


.0006 


— .0021 






— .0028 


—.0027 


— .0027 








.0001 


.0000 










— .0001 













A23/0 

A> 
A^'^yo 

yx=yo-\-x^yo-\ — , — '-^^yo+^ 1^^ '-^^yo 

|2 |3 



x{x — i) (a; — 2) (x — 3) 



A*yo+etc. 



= -.i826+|(-.2637)- 



K-f) 



(.0061)- 



K-f)(-l) 



(-.0028) 



l(-f)(-i)(-¥) 
24 



(.0001) 



= -•3417- 



3. The following example is taken from Boole's Finite Differ- 
ences. Given log 3. 14 = .4969296, log 3.15 = .4983106, log 3.16 = 
.4996871, log 3. 17 = .5010593; required an approximate value of 
log 3.14159- 

yo yi y-i yz 



Ayo- 

A2jo. 

A^yo. 



.4969296 .4983106 .4996871 .5010593 
.0013810 .0013765 .0013722 
— . 0000045 ~ ■ 0000043 
. 0000002 

Here the value of x in the formula is equal to 0.159. 
3'x = .4969296+(-i59)(-ooi38io)+-!^i^2ZkI59 — L!(_.ooooo45) 



•i59(-i59-i)(-i59-2) 



(.0000002) 



= .4971495. 
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This is correct to the last decimal place. If only two terms 
had been used in the right-hand member of the formula, which 
is equivalent to the rule of proportional parts, there would 
have been an error of 3 in the last decimal place. The rapid 
decrease in the value of the differences enables us to judge 
quite well of the accuracy of the results. The above formula 
can be applied only when the values of x form an arithmetical 
series. 

In case the series of values given are not equidistant, that is, 
the values of the independent variable do not form an arithmetical 
series, it becomes necessary to apply another formula. 

Let ya, ji, yc, ya, ■ ■ • yt. be the given values corresponding 
to a, b, c, d, . . . k respectively as values of x. It is required 
to find an approximate expression for yx, an unknown term 
corresponding to a value of x between x = a and x = k. 

Since there are n conditions to be satisfied the expression 
which is to represent all of the values must contain n constants. 
Assimie as the general expression 



y:,=A+Bx-\-Cx^+Do(?+ . . . +N3(f'-\ 

Geometrically this is equivalent to drawing through the n 
points represented by the n sets of corresponding values a 
parabola of degree w— i. 

Substituting the sets of values given by the data in the 
equation above n equations are obtained from which to determine 
the values oi A, B, C, etc., 

ya=A+Ba+Ca^+Da3+ . . . Na"-^; 
yi=A+Bb+Cb^+Db^+ . . . Nb"'^; 



yi=A+Bk+Ck^+DP+ . . . Nk"-\ 

But the solution of these equations would require a great 
deal of work which can be avoided by using another but equiva- 
lent form of equation. 
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Let y;,=A{x-b)(x-c)(x-d) . . . (x-k) 

+B{x-a){x-c){x-d) . . . (x-k) 

+C{x-a)(x-b)ix-d) . . . (x-k) 

+D(x-a)(x-b)(x-c) . . . (x-k) 
+ etc. to n terms. 

Each one of the n terms on the right-hand side of the equation 
lacks one of the factors x—a, x—b, x—c, x—d, . . . x-k, 
and each is affected with an arbitrary constant. The expression 
on the right-hand side of the equation is a rational integral 
function of x. 

Letting x = a gives 



and 



ya=A(a — b)(a—c)(a—d) 



a—k. 



A=- 



ya 



(a — b)(a — c)(a—d} . . . a — k' 
Letting x = b gives 
B=- 



yo 



(b-a)(b-c)(b-d) . . . (b-k)' 

Proceeding in the same way we obtain values for all of the 
constants and, finally, 

(x — b)(x — c)(x — d) . . . (x — k) 



yx=ya- 



+yi 



+y 



(a — b)(a — c)(a — d) 
(x — a)(x — c)(x—d) 



(b-a)(h-c)(b-d) 
(x — a) (x — b) (x — d) 



+yd 



(c — a)(c —b)(c — d) 
(x — a)(x — b)(x — c) 



(d-a)(d-b)(d-c) 



+yr 



(x — a)(x — b)(x — c) 



(k-a)(k-b)(k-c) . 



. (a-k) 
. (x-k) 



. (b-k) 
. (x-k) 



. (c-k) 
. (x-k) 



. (d-k) 



(2) 
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This is called Lagrange's theorem for interpolation. 

I. Apply formula (2) to the data given under formula XIX 
for finding the value of y corresponding to a; = 0.425. Select 
two values on either side of the value required, 





a^-Z^, > = 


•69S> 






&=,40, > = 


•73°; 






c = .4S, yc = 


•757, 






i = -$o, ya = 


.780. 




a; in 


the formula must be taken as c 


>-S- 






-C«(^^^^^^ 


,...,a(-*)(-t) 




, (..^^^-^ 


+ (".-780 


Mm-V) 



'(2)(l)(-l) 



(3)(2)(i) 



= .744. 



2. Required an approximate value of log 212 from the fol- 
lowing data: 

log 210 = 2.3222193, 

log 211 = 2.3242825, 

log 213 = 2.3283796, 

log 214 = 2.3304138. 



log 212 = (2.3222193) 
+ (2.3283796) 



(i)(-i)(-2) 



-(2.3242825) 



(2)(-l)(-2) 



(-i)(-3)(-4) ' " " "'(i)(-2)(-3) 

(^)(i)(-^) +(,3304X38)i^Kl)(^) 
+ U.3304i3«J (^)(^)(^) 



(3)(2)(-i) 
= 2.326359. 

This is correct to the last figure. 

, In case the values given are periodic it is better to use a 
formula involving circular functions. In Chapter V the approxi- 
mate values of the constants in formula XX were derived. This 
formula could be used as an interpolation formula. But on 
account of the work involved in determining the constants it is 
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much more convenient to use an equivalent one which does not 
necessitate the determination of constants.* The equivalent 
formula given by Gauss is 



yx=ya 



sin |(x— J) sin \{x—c). 



sin|(a;— ^) 



sin §(a — i) sin |(a — c) . . . sin|(a — yfe) 

sin ^(:k— a) sin I (a;— c) . . . sin|(x— ^) 

sin 1(6 — a) sin |(&— c) . . . sin|(&— ^) 

sin J(a; — a) sin |(x— S) . . . sin^(x— ^) 



sin \{c — (i) sin \{<^ — h) 
+ etc 



sin \{(^—^ 



(3) 



It is evident that the value of ja is obtained from this formula 
by putting x = a. The value of y„ is obtained by putting x = h, 
and jc by putting x = c. 

The proof that (3) is equivalent to XX need not be given 
here. 

Let it be required to find an approximate value of y cor- 
responding to a; = 42° from the values given. 



From (3) 
3» = (io.i) - 



X 


y 


30° 


10. 1 


40" 


9.8 


50" 


8.5 



sin I ° sin ( —4°) 



.sin (—5°) sin( — 10°) 



-(9-8) 



sin 6° sin (—4°) 
sin 5° sin (-5°) 



+ (8.S) 4 



sin 6° sin i' 



sin 10 sin 5 



, . (.0175) (.0698) , . (.1 045) (.0698) 
■(^°"\.d872)(.x736) + ^^-'^~(.o872)^ 

(.1045) (.0175) 



+(8.5) 



= 9.618. 



(.i736)(.o872) 



* Trigometrische Interpolation, Encyklopadie der Mathematischen 
Wissenchaften, Vol. II, pt. I, pp. 642-693. 
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A better result would have been obtained by using four sets 
of values. 

Differentiation of Tabulated Functions 

It is frequently desirable to obtain the first and second 
derivatives of a tabulated function to a closer approximation 
than graphical methods will yield. For that purpose we will 
derive differentiation formulas from (i) and (2). From 

, ^ ,x(x-l)^„ , x(x-l)(x — 2)^^ 

jx = yQ-\-x^yQ +-^-| — -t^yo-\- , -l^ya 

^ x{x-t){x-2){x-i) ^, _l_ 

k 

, x{x-\){x-2) . . . (x-n+2)^ „ _ 1 , . 



n — 1 



By differentiating it follows that 



y'. = Ayo +^^^A2yo+2^^-^^^±^A3yo 

2 Is 

^ 4X^-.2X^ + 22X-6 ^,^^^ ^^^ 

|4 
Differentiating again 
yx" = A''yo+{x-i)A^yo+(^x^-x+ii)A'yo+ (S) 

As an illustration let it be required to find the first and 
second derivatives of the function given in the table below and 
determine whether the series of observations is periodic* 

The consecutive daily observations of a function being 
0.099833, 0.208460, 0.314566, 0.416871, 0.514136, 0.605186, 
0.688921, 0.764329, show that the function is periodic and deter- 
mine its period. 

* Interpolation and Numerical Integration, by David Gibb. 
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From the given observations the following table may be 



written: 

% y =/(«) 

1 0.099833 

2 0.208460 
0.314566 
O.416871 
0.514136 
0.605186 
0.688921 



(- 
(- 
(- 



8 0.764329 

From (4) 

/i= .108627 
.001260 

. 109887 
-.000437 

• 109450 

y'z = . 102305 
.002520 



. 104825 
.000411 

. 104414 



0.108627 
0.106106 
0.102305 
)o.o97265 
)o.o9io5o 
)o.o83735 
)o.o754o8 



— .000427 
— .000010 

-.000437 



— .000392 
— .000019 

— .000411 



.002521 
.003801 
.005040 
.006215 

•0073 IS 
.008327 



A' 

— .001280 
(+) -.001239 

(+)-.ooii75 
(+) — .001100 
(+) — .001012 



A< 



.000041 
.000064 
.000075 
.000088 



y2-- 



.106106 
.001900 

. 108006 
— .000429 

.107577 

y'i= .097265 
.003108 



• 100373 
— .000389 

■099984 



.000413 
.000016 

.000429 



.000367 
.000022 

. 000389 



For the remaining first derivatives the order must be reversed 
and the resulting sign changed. 

3;'5= — .097265 .002520 y6= — .091050 
.000413 



0000 10 



.000016 



.003108 
.000392 



097275 
002933 

094342 



.002933 



.091066 
.003500 

.087566 



.003500 



no 
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/7=-. 083 73 5 

— .000019 



083754 
004025 



— .006279 
. 000038 

— .006241 

y"7=-.oo73is 
— .001100 



— .008415 
. 000069 

— .008346 



.003658 ^8= -.075408 
.000367 —.000022 



.004025 



.079729 

From (5) 






y"\ = — .002521 


.001280 


y"2 


.001318 


.000038 
.001318 




— .001203 




/';.,= -.005040 


.001175 


y'4 


.001244 


. 000069 
.001244 




-.003796 




y"h = — . 005040 




y"e 


— .001239 







075430 

004501 



y 8= 



070929 



003801 
001298 



002503 

006215 
001181 



005034 

006215 
00117s 



007390 
000059 



007331 

008327 
001012 



009339 
000081 



009258 



.004164 
■000337 

.004501 



.001239 
. 000059 

.001298 

.001100 
.000081 

.001181 



X 


y 


y' 


y" 


y 
— .0121 


I 


•099833 


. 109450 


— .001203 


2 


. 208460 


.107577 


-.002503 


— .0120 


3 


.314566 


. 104414 


--.003796 


— .0121 


4 


.416871 


.099984 


-.005034 


— .0121 
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X 


y 


y' 


y" 


yl 


5 
6 

7 
8 


•514136 
.605186 
.688921 
.764329 


.094342 

.087566 - 
.079729 
.070929 — 


.006241 
.007331 
. 008346 
.009258 


— .0121 
— .0121 
— .0121 
— .0121 



Since ^^ IS very nearly constant and equal to -.0121, the 
corresponding differential equation is 

y"+.oi2iy = o, 
whose solution is 

y = A coso.iia;+5 sino.iia;. 

/ 

This shows that ;y is a period function of x, and its period is 

2ir , 

■ — , or 57.12 days. 

Convenient formulas for the first and second derivatives may 
also be obtained by differentiating Lagrange's formula for inter- 
polation. Using five points the ^formula is 

_ {x—b)(x—c){x—d)(x—e) (x—a){x—c)(x—d)(x—e) 
•^'~^'' {a-b){a-c){a-d)(a-ey^'' {b-a)(b-c){b-d){b-e) 

. (x — a){x — b)(x — d)(x — e) (x — a)(x — b){x — c)(x — e) 
■""^^ {c-a)ic-b)(c-d){c-e) '^^'^ {d-a){d-b){d-c){d-e) 

{x-a){x-b){x-c){x-d) , . 

"^■^' {e-a){e-b){e-c){e-d) ^ ' 

Selecting the points at equal intervals and letting 

e—d=d—c = c — b = b—a=h, 

and differentiation 

y'a=^[-25ya+4Sy!,-3^yc+T-6yd-3ye], 

12 n 

y'b = — A- 3ya-ioyb+i?>yo- 6ya+ ye], 

12/2 
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y'c=—T[ >- 83/6 

I2« 



+ Sya-ye], 



y'a^^[- ya+ 6yt,-i8yo+ioya+3ye], 

y'e=-^[ sya-^^yt+s^yc-A&ya+^syel 

Differentiating again 

y"a = -^bsyo - 10476+ 1 uyc - 567*+ 1 lyJ, 

y'i,=-^[ii3'a- 20^6+ 6yc+ 4yd-ye], 
i2h'' 

y"c=-^{-ya-{- i6ys- 303/.+ i()yd-ye], 

y'd=-^{-ya+ 43'!-+ ^y'>- 2oya-\-iiye], 
y'e = -^[i 13'a - 563'*+ 1 143'c- i043'<i+3S)'«]- 

I2«'' 

The results of applying these formulas to the function given 
are expressed in the table below. 



X 


>> 


y 


y" 


I 


■099833 


. 10945 I 


001203 


2 


. 208460 


• 107583 


002524 


3 


.314566 


.104415 


003804 


4 


.416871 


. 099986 — 


005045 


5 


•514136 


•094347 


006221 


6 


.605186 


.087568 - 


007322 


7 


.688921 


■079733 


008334 


8 


.764329 

_1j__ J*. ■„!__ 


.070929 — 

n 'ii^ ii 


009258 

1 1 



These results agree fairly well with those previously obtained. 
It is probable that the formulas derived from the interpolation 
formula give the most satisfactory results. 
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As another application let us find the maximum or minimum 
value of a function having given three values near the critical 
point. 

Let ja, yt, and jc be three values of a function of x nea-r 
its maximum or minimum corresponding to the values of x, 
a, b, and c respectively. 
From (2) 

_ {x — b)(x — c) (x—a)(x—c) (x—a)(x—b) 
y^-y" (^a-b){a-cy^\b-a){b-cy^' {c-a){c-by 

Equating to zero the first derivative with respect to x 

2x—b—c , 2x—a—c . 2x — a — b 



y^=ya 



-y<, 



-yc- 



=0; 



{a-b){a-c) ' ''" {b-a){b-c) ' ^' {c-a){c-b) 

^ _ y.(b^ - c^) +yt{c^ - a^) +yc(a^ - b^) (6) 

2\ya{b-c)+yb{c-a)+yc{a-b)] 

This is equivalent to drawing the parabola 

y=A-\-Bx-\-Cx'^ 

through the three points and determining its maximum or 
minimum. 



From the table of values 






X 


^ 


; 


6.0 


10 


OS 


6.5 


10 


14 


7.0 


10 


10 



the abscissa of the maximum point is found from (6). 

_ (io.o5)(-6.7.'i) + (io-i4)(i.s) + (io-io)(-6-25) _g ^5 
"" 2[(io.o5)(-.5) + (io.i4)(i) + (io.io)(-.5) 



y = 10.1424. 



CHAPTER VIII 

NUMERICAL INTEGRATION 

Areas 

An area bounded by the curve, y=f(x), the axis of x, and 
two given ordinates is represented by the definite integral 



-£■ 



ydx, 

where the ordinates are taken at x = a and x=n. It may be 
said that the definite integral represents the area under the 
curve, or that the area under the curve represents the value of 
the definite integral. 

If a function is given by its graph, it is possible, by means 
of the planimeter, to find roughly the area bounded by the curve, 
two given ordinates and the a; — axis, or, what amounts to the 
same thing, the area enclosed by a curve. This method is used 
in finding the area of the indicator diagrams of steam, gas or 
oil engines, and various other diagrams. The approximations 
in these cases are close enough to satisfy the requirements. 

If, however, considerable accuracy is sought, or whenever 
the function is defined by a table of numerical values another 
method must be employed. 

Mechanical Quadrature or Numerical Integration is the method 
of evaluating the definite integral of a function when the func- 
tion is given by a series of numerical values. Even when the 
function is defined by an analytical expression but which can- 
not be integrated in terms of known functions by the method 
of the integral calculus, numerical integration must be resorted 
to for its evaluation. 

The formulas employed in numerical integration are derived 
from those established for interpolation. 

114 
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In interpolation it was found that the order of differences 
•which mUst be taken into account depends upon the rapidity 
with which the differences decrease as the order increases. 
This is also true of numerical integration. It is the same as 
saying that if the series employed does not converge the process 
will give unsatisfactory results. An illustration will be given 
later. 

Formulas for numerical integration will be derived from (i) 
of Chapter VII. 

In this formula it was assumed that the ordinates are given 
at equal intervals. 

yx=yo+xAyo+-^-. -A^yo+ — /^ -A^yo 

11 \a 

, x{x-i)(x-2)(x-7,) ^^^^ , x(x-i)(x-2){x-s)(x-4) ^^^^ 
H i A^yo H i A°yo 

li Is 

^ xix-i)(x-2)ix-s)(x-4)(x-s) ^^ ^ ^ ^ . . (i^ 

Integrating the right-hand member, 
I yi.dx=yo \,dx+Ayo\ xdx-\ — p- | x{x — i)dx 

Jo Jo \2_J 

-\ — 1^ I x(x — i){x — 2)dx 
13 Jo 

+^ {"xix-i){x-2)ix- 2)dx 
W Jo 

+^rx{x-i){x-2){x-3)(x-4)dx 
+^rx{x-i)ix-2){x-3){x-4Kx-5)dx+ . . . 



+ ---'— +^^-3^ 



/T 
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\6 4 3 /Is 

\7 2 ' 4 3 / P 

The data given in any particular problem will enable us 
to compute the successive differences of yo up to A"yo. On the 
assumption that all succeeding differences are so small as to 
be negligible the above formula gives an approximate value 
of the integral. It is only necessary to assign particular values 
to n. 

Let w = 2, then 

I y^x = 2yo+2Ayo+iA^yo, 

Ayo=yi-yo, 

A^yo=Ayi-Ayo=y2-yi -yi+yo, 
=y2-2yi+yo. 

Substituting these values in the above integral it becomes 
yxdx = 2yo+2yi - 2yo+iy2—%yi+iyo, 



j: 



_ yo+4yi+y2 

3 



This is equivalent to assuming that the curve coincides with 
a parabola of the second degree. 

If the common distance between the ordinates is h, the 
value becomes 

ydx=ih(yo+4yi+y2) (7) 



Jr2i 




ydx=3yo+iAyo+iA^yo+^A^yo, 

Ayo=yi-yo, 
A^yo =Ayi-Ayo=y2- 2yi +yo, 



i 



NUMERICAL INTEGRATION 117 

A^yo = A^yi - A^yo = ^y•2. - ^y\ - Ayi + Ayo 
=y3-33'2+ 3^1 -Jo- 
Substituting these values in the equations, 

I y^a; = 3yo+l:yi-l5'o+f3'2-f3'i+f3'o+l3'3-l3'2+l3'i-b'o, 

= l3'o+l3'i+l0'2+l3'3, 
=f(yo+33'i+33'2+3'3)- 

If the common distance between the ordinates is h the 
formula becomes 

= tA(3'o+33'i+33'2+J'3) (8) 







This is equivalent to assuming that the curve coincides with 
a parabola of the third degree. 

If there are five equidistant ordinates, h representing the 
distance between successive ordinates 



i 



45 



If the area is divided into six parts bounded by seven equi- 
distant ordinates the integral becomes 

I 3;Ja; = 6yo + i8A3/o+27A2);o+24A33;o+ WA^jo 
+ffA53'o+AVA>. 

Sine 3 the last coefficient, AV, differs but slightly from -h 
and by the assumption that A^j^o is small the error will be slight 
if the last coefficient is replaced by -i^. 

Doing this and replacing 

Lya by y\ — yo, 
A^j^oby )'2 — a^/i+yo, 
A^j^o by y3-3y2+3:)'i-3'o> 
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A^yo by yi—4y3+(>y2 —4yi+yo, 
A^yo by ys - 5^4 +io3'3- 103/2 +5^1 —Jo, 
A^yo by y^ — 6^5+1 5^4 — 20^3 + 1 5^2 — 63^1 +^0, 
gives the formula 



Jr6h 
I' 



ydx = ^h\yo+y2+yi+y6+5(yi+y5)+6y3\- • i^°) 




The application of these 
formulas is illustrated by 
finding the area in Fig. 27. 



Fig. 27. 

By (7) 
A=ih(yo+4yi+2y2+4y3+2yi+4y5+2y6+4y7+2ys+4yg 
+2ym+4yn+yi2). 

By (8) 
A=U(yo+3yi+3y2+2y3+sy4+sy5+2y6+3y7+sys+2y9 

+3yio+3yii+yi2)- 

By (9) 

A=':hh[i4(yo+2y^+2ys+yi2)+64(yi+y3+y5+y7+y9+yn) 

+24(y2+y6+yio)]- 
By (10) 

A=^h\yo+y2+yi+2y6+y8+yio+yi2+5(yi+y5+y7+yii) 

+(>(y3+y9)]- 

I. A rough comparison of the approximations by the use 
of these formulas will be obtained by finding the' value of 

— . The value of this definite integral is log 13 = 2.565. It 
1 X 

is also equal to the area under the curve 

I 

y=- 

X 
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from x = i to a;-=i3. Dividing the area up into 12 strips of 
unit width by 13 ordinates the corresponding values of x and 
y are 



X 


I 


2 


3 


4 


S 


6 


7 


8 


9 


10 


II 


12 


13 


y 


I 


1 

"2" 


i 


i 


i 


i 


+ 


i 


i 


1 


1 
1 1 


^ 


tV 



By (7) 

^=Mi+2+f+i+l+l+f+Hf+t+A+i+TV] 

= 2.578, error .5%; 

By (8), ^ = 2.585, error .8%; 

By (9), 4 = 2.573, error .3%; 

By (10), A = 2.572, error .3%. 

2. The accuracy of the approximation is much increased 
by taking the ordinates nearer together, as is shown by the 
following evaluation of 

X' 

The value of this integral is equal to the area under the 
curve 



dx 
i+x' 



I 

y=- 



i+x' 

from x=o ta x = i. Dividing the area into twelve parts by 

— — is found to be 

By (7), 0.69314866, error 0.00000148; 

By (8), 0.69315046, error 0.00000328; 

By (9), 0.69314725, error 0.00000007; 

By (10), 0.69314722, error 0.00000004. 

The correct value is, of course, log, 2, which is 0.69314718. 
Formulas (7) and (8) are Simpson's Rules, (10) is Weddle's 
Rule. 
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3. Apply the above formulas to the area of that part of the 
semi-ellipse included between the two perpendiculars erected at 
the middle points of the semi-major axes. Let this area be 
divided into twelve parts by equidistant ordinates. 

Since the equation of the elUpse is 

these ordinates are 

2V3 &, i^V^g h, |V8 I, \V7i b, iV3S b, tJ^Vhs 6, b, 
T^V^ b, WYs b, |V7^ b, \Vl b, T^jV^ b, IV3 b. 

By (7), ^=o.9566o99a&; 

By (8), ^=o.9566o8oa&; 

By (9), ^=o.9566ii4a6; 

By (10), A =0.956611406. 

The correct value to seven 
places is 0.956611506. 

In the appHcation of these 
formulas it is highly desirable 
to avoid large differences among 
the ordinates. For that reason 
the formulas do not give so 
good results when applied to 
the quadrant of the ellipse. 

4. The area under the curve, 
Fig. 28, determined by the 
following sets of values: 




A .6 .8 1.0 1.2 a 

Fig. 28. 



.2 



■4 



.6 



i.o 



1.2 



2.2 



2.7 2.6 



2-3 



is by (7) 

^=i--Ki-o+6.o-l-44+io-8+S-2+9.2-h2.i) = 2.s8, 
and by (8), 

^ =1-1(1.0-1-4.5-1-6.6-1-5.4 -f7.8H-6.9+2.i) = 2.5725. 



2.1 
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This area is represented by the definite integral ( ydx 



I" 



The area found is therefore the approximate value of this 
integral 

5. Find the area under the curve determined by the points 

x|i 1-5 1-9 2.3 2.8 3-2 3.6 4.0 4-6 4.8 5.0 
y\o .40 1.08 1.82 2.06 2.20 2.30 2.25 2.00 1.80 1.5 

The points located by the above sets of values are plotted 
in Fig. 29 and a smooth curve drawn through them. The area 
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Fig. 29. 

is divided into strips each having a width of .4. Rectangles 
are formed with the sanie area as the corresponding strips. 
The eye is a very good judge of the position of the upper bound- 
ary of each rectangle. Adding the lengths of these rectangles 
and multiplying the sum by .4 the area is found to be 6.644. 
By Simpson's Rule, formula (7), are found 



for h = .2, 
h = A, 



4=6.639, 
4=6.645. 



The graphical determination of areas can be made to yield 
a close approximation by taking narrow strips, and where the 
points are given at irregular intervals the area can be found 
more rapidly than by the application of Simpson's Rules. 
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6. A gas expands from volume 2 to volume 10, so that its 
pressure p and volume v satisfy the equation pv = ioo. Find 
the average pressure between v = 2 and v = 10. 

The average pressure is equal to the work done divided by 

8. The work is equal to the area under the curve p = — from 

' V 

v = 2 to D = 10, which is 

ri«ioo, . r, V" A 

I — dv = 100 log V = 160.944. 
J2 i> I h 

That this area represents the .work done in expanding the 
volume from 2 to 10 becomes evident in the following way. 
Let s represent the surface inclosing the gas, ps will then be 
the total pressure on that surface. The element of work will 
then be 

dW=psdn, 

when dn lepresents the element along the normal. 



= j psdn. 



W 

But 

sdn=dv, 
and 

TF= ( pdv. 

This is the equation above. The average pressure over the 
\. change of volume from 2 to 10 is 

160.944 -i-8 = 20.118. 

7. Find the mean value of sin^ x from x=o to a: = 2x. Plot 
the curve y = sm^ x by the following values of x and y: 



X 





T 
12 


6 


TT 
4 


T 

3 


51 
12 


TT 
2 


y 





.0670 


.2500 


.5000 


.7500 


•9330 


I. 0000 



X 


11 
12 


27r 
3 


31 
4 


St 
6 


I IX 
12 


y 


• 9330 


.7500 


.5000 


.2500 


.0670 
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I37r 
12 



6 



4 



4Tr 
3 



I77r 

12 



31 

2 



y 


o 


.0670 


.2500 


.5000 


.7500 


■9330 


I . 0000 


X 


IQTr 
12 


Si 
3 


Ttt 
4 


IITT 
6 


2_37r 
12 


2ir 




y 


•933° 


• 7500 


.5000 


.2500 


.0670 








Applying Simpson's Rule, formula (7), the area is found to 
be TV. The mean value is the area divided by 27r or .5. 

8. A body weighing 100 lb. moves along a straight line 
without rotating, so that its velocity v at time t is given by the 
following table: 

^ sec 



V ft./sec 



1.47 1.58 1.67 1.76 1.86 

Find the mean value of its kinetic energy from t = i to t = g. 



t 


I 


3 


5 


7 


9 


li" 


2 . 1609 


2.4964 


2 . 7889 


3.0976 


3-4596 


Kinetic energy . 


3-355 


3.876 


4331 


4.810 


5-372 



Plotting kinetic energy to t, the area under the curve is 
34.755. This divided by 8 gives the mean kinetic energy as 

4-3S7- 

Volumes 
Fig. 30 explains the ap- 
phcation of the formulas 
to the problem of finding 
the approximate volume of 
an irregular figure. Tiie area 
of the sections at right angles 
to the axis of x are : 

Ai=\k{yi+Ay5-\-y'd, 
A2 = \k{y&-[-^y^-\-ys), 
Az=\k{y2-^Ay7+y^)- Fig. 30. 
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If the areas of these sections be looked upon as ordinates, 
h being the distance between two adjacent ones, it is evident 
that the volume may be represented by the area under the 
curve drawn through the extremities of these ordinates. 

V=\KAi+aA2+Az) 

Substituting the values of Ai, A2, and A3 in this equation, 
the volume becomes 

v=\h[^k{yi+4y5+yd+U{y6+4y9+y8)+U(y2+4y7+y3)] 
=ihk[yi+y2+y3+yi+4(y5+y6+y7+y8)+i6yQ] 

In order to apply formulas (8), (9) and (10), the solid would 
have to be divided differently, but the method of application 
is at once evident from the above and needs no further discussion. 

1. The following are values of the area in square feet of the 
cross-section of a railway cutting taken at intervals of 6 ft. 
How many cubic feet of earth must be removed in making the 
cutting between the two end sections given? 

91, 95, 100, 102, 98, 90, 79. 

These cross-section areas were obtained by the application 
of Simpson's Rules. 
By (7), 

F = f -6(91 -1-380 -I- 200 -1-408 -F 196 4-360-1-79) =3428; 

By (8), 

F = |-6(9i-F28s-F30o-l-204-f 294 -f-27o-F79) =3426.8. 

2. .4 is the area of the surface of the water in a reservoir 
when full to a depth h. 



hit.... 


30 


25 


20 


15 


10 


5 





A sq.ft. . 


26,700 


22,400 


19,000 


16,500 


14,000 


10,000 


S,ooo 
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Find (a) the volume of water in the reservoir, (b) the work 
done in pumping water out of the reservoir to a height of loo ft. 
above the bottom until the remaining water has a depth of 
loft. 

^^=#(26,700+89,600+38,000+66,000+28,000+40,000+5,000) 
=488,833 cu. ft. 

rzo 
Work=w I A{ioo-h)dh, where w= weight of i cu.ft. of water 

= 62.3 lb. The value of this integral will be approximately the 
area under the curve determined by the points 



h 


30 


25 


20 


15 10 


A{ioo — h). 


1,869,000 


1,680,000 


1,520,000 


1,402,500 1,260,000 



multiplied by 62.3. 
This area is equal to 

1(1,869,000+6,920,000+3,040,000+5,610,000+1,260,000) 
= 31,165,000. 

Multiplying this by 62.3 gives the work equal to 1,941,579,500 
ft.-lb. 

3. When the curve in Fig. 29 revolves about the ic-axis, 
find the volume generated. 

The areas of the cross-sections corresponding to the given 
values of x are given in the following table: 

o .2 .4 .6 .8 i.o 1.2 



A.... IT 2.257r 4.84X 7.2977 6.7677 5.2977 4.4177 

By (7) 7 = 5.86277=18.416. 

By (8) 7 = 5.80377=18.231. 

4. When the curve in Fig. 30 revolves about the a;-axis, 
find the volume generated from x = i tox=4.2. From the curve 
the following sets of values are obtained : 
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X i.o 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 



. II 



•29 -53 



J_ 

y^ o .012 .084 .281 
X 2.8 3.0 3.2 3 



87 1.37 I. 71 1.90 



2.01 



757 1.877 2.924 3.610 4.040 
4 3-6 3-8 4-0 4-2 



y 



2 .06 2.12 2.2 



27 2.30 2.28 2.25 



2.20 



y-' 4.244 4.494 4.84 5.153 5.290 5.198 5.062 4.84 
The volume is by (7) 

T-i-i(i49.oo4) =31.2 cu. units. 

Centroids 
Let the coordinates of the centroid of an area be represented 
by X and y. Then from the calculus 






xydx 



C" 
I C 

- \ ^y^dx 
I ydx 



The integral in the numerator of the value of i may be 
represented by the area boimded by the curve Y=xy, the x-axis 
and the two ordinates x = a and x = b. The original area is 
bounded by the curve whose ordinates are represented by y, 
the X-axis and the two ordinates x = a and x = b. The integral 
in the numerator of the value of y may be represented by the 
area bounded by the curve Y=y^, the x-axis and the two 
ordinates x = a and x = b. 

For a volume generated by revolving a given area about the 
X-axis 






y^xdx 



y^dx 
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Axdx 



Adx 



The process of finding the coordinates of the centroid of the 
area in Fig. 28 is shown in the table: 



X 





.2 


■4 


.6 


.8 


1.0 


1.2 


y 


I.O 


i-S 


2.2 


2.7 


2.6 


2-3 


2.1 


xy 


0.00 


0.30 


0.88 


1.62 


2.08 


2.30 


2.52 


J,2 


1. 00 


2.25 


4.84 


7.29 


6.76 


S-29 


4.41 


y'x 


0.000 


0.450 


1.936 


4-374 


5.408 


S.290 


5.292 



The area under the curve Y=xy is 
tf[o.oo+i.2oH-i.76+6.48+4.i6+9.2o+2.52] = 1.688; 
- 1.688 



2.58 



= .654. 



The area under the curve Y=^y^ is 
^[1.00+9.00+9.68+29.16+13.52+21.16+4.41] = 2.931 

^ = -^=1.136. 
2.58 

As was pointed out before, large changes in the ordinates 
must be avoided. 

For the volume generated by revolving the area about 
the a;-axis 

-_ 7rA[o-ooo+i-8oo+3.872 + i7.496+io.8i6+2i.i6o+5.292] 
TtV[i-oo+9.oo+9.68+29.i6+i3.52+2i.i6+4.4i] 

^60436^68 
87-93 
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Moments of Inerila. 

The expression for the moment of inertia of an area about 
the 3/-axis is 



About the a;-axis 



C" 
iy= I x^ydx. 

c= I xy^dy. 



When the equation of the curve is known these integrals 
can be calculated at once, but when this is not the case approxi- 
mate methods must be resorted to. 

I. The process of finding the approximate values of these 
integrals is shown in the table below. The values of x and y 
are taken from Fig. 28. 



X 





.2 


.4 


.6 


.8 


1.0 


1.2 


y 


I.O 


i-S 


2.2 


2.7 


2.6 


2-3 


2. 1 


x^y 


0.000 


0.060 


0-33S 


0.972 


1.664 


2.300 


3-024 


h' 


0-333 


I.I2S 


3-549 


6.561 


5-859 


4.056 


3.087 



If the values of x^y be plotted to x we will have a curve 
under which the area represents the moment of inertia of the 
area in Fig. 28 about the ji-axis. 

A = i-357- 

Dividing this by the area found before, there results for the 
radius of gyration about the y-axis 

Plotting ^'^ to X and finding the area imder the curve so 
determined 

7:r =4-6136, 
and 

i?x2 = 1.788. 
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2. The form of a quarter section of a hollow pillar, Fig. 31, 
is given by the following table. Find the moment of inertia 
of the section about the axes of x and y. 



.5 


i 
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Fig. 31. 



X 


Y 


xW 


y 


X 


y^X 


.00 


• 050 


. 00000 


.00 


.100 


. 00000 


•OS 


•OSS 


.00014 


•OS 


.108 


.00027 


.10 . 


.068 


.00068 


.10 


.116 


.00116 


• IS 


.078 


.00175 


•IS 


.120 


.00270 


. 20 


.096 


•00384 


. 20 


•125 


. 00500 


•25 


.116 


.0072s 


•25 


.130 


.00812 


• 30 


.148 


•01332 


• 30 


• 133 


.01197 


■iS 


.200 


.02450 


•3S 


.140 


•0171S 


.40 


• 300 


. 04800 


.40 


•ISO 


.02400 


■AS 


•21S 


.043S4 


•45 


•215 


•04354 


■so 


.000 


. 00000 


• SO 


.000 


. 00000 



In the above table X stands for the width of the area parallel 
to the X-axis and Y for the width parallel to the y-axis. The 
area is 0.066. 

The moment of inertia about the y-axx?, is 



X 



a:2p'c?x = . 00736; 

.066 
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The moment of inertia about the x-axis is 

Jy^Xdy = .oo6ig; 


„ „ .oo6ig 

where R stands for the racMus of gyration. 

The values of the above integrals were computed by for- 
mula (7). 
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If a chart could be constructed with all the different forms 
of curves together with their equations which may arise in 
representing different sets of data it would be a comparatively- 
simple matter to select from the curves so constructed the one 
best suited for any particular set. Useful as such a chart would 
be its construction is clearly out of the question. The most 
that can be done of such a nature is to draw a number of curves 
represented by each one of the simpler equations. 

A word of caution is, however, riecessary here. A particular 
curve may seem to the eye to be the one best suited for a given 
set of data, and yet, when the test is applied, it may be found 
to be a very poor fit. It is of some aid, nevertheless, to have 
before the eye a few of the curves represented by a given formula. 

The purpose of the following figures is to illustrate the 
changes in the form of curves produced by slight changes in 
the constants. Figs. I, II, III, and IV show changes produced 
by the addition of terms, Figs. V to XIX changes in form 
produced by changes in the values of the constants, and Fig. 
XX the changes in form brought about by varying both the 
values of the constants and the munber of terms. 

A discussion of all the figures is unnecessary. A few words 
in regard to one will sufi&ce. Formula XIV, for example, 
y = a+bx', an equation which can be made to express fairly 
well the quantity of water flowing in many streams if x 
stands for mean depth and y for the discharge per second, 
represents a family of triply infinite number of curves. Fixing 
the values of b and c and varying the value of a does not 
change the form of the curve, but only moves it up or down 
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along the y-axis. Keeping the values of a and b constant and 
varying the value of c, the formula will represent an infinite 
number of curves all cutting the jz-axis in the same point. In 
the same way, keeping the value of a and c constant and vary- 
ing the value of b, an infinite number of curves is obtained, 
all of which cut the jz-axis in a fixed point. In Fig. XIV the 
quantity a is constant and equal to unity, while b and c vary. 

To one trained in the theory of curves the illustrations are, 
of course, of no essential value, but to one not so trained they 
may be of considerable help. 

The text should be consulted in connection with the curves 
in any figure. The figures are designated to correspond to the 
formulas discussed in the first five chapters. 
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(1) y=\-.\x 

(2) !/=l-.la;-f.01i2 

(3) i/=l-.l2:+.01i2-,001x= 

(4) !/=l-.l3:+.01i2_.001i'-|- 

.OOOll* 

(5) i/=l-.ls-l-.01i=-.001x'-|- 

.OOOli^-.OOOOla;* 

(6) ;/=-l-.l^-|-.01a;2-.001x'+ 

.OOOlrc^-.OOOOli^-t-.OOOOOlx^ 

See formula I, page 13 
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(1) 8(=1-1A 

(2) y^l-l/x+l/x" 

(3) !/=l-l/i+l/x2-l/a:3 

(4) !/=l-l/i-M/i2-l/i3+l/a:4 

(5) y=l-l/z+l/x'~l/x'+ 

l/x*-l/x^ 

(6) !/=l-l/i-l-l/j2_i/^3+ 

l/x*~l/x^+l/x^ 
See formula II, page 22 



Fig. II. 




(1) 
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.Ix 




(2) 
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=1- 


.li-l-.01a;2 




■,3) 


1 
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=1- 


.li-F.OU2_ 


.OOli' 


(4) 
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=1- 


.li-|-.01i2- 


.00U3_,_ 




.000 Ix* 





{£) -=l-.li+.01x2_.ooix3+ 
.OOOli^-.OOOOlx^ 

(6) - = l-.li+.01x2-. 001x5+ , 
V 
.OOOlx^-.OOOOlx^+.OOOOOla;* 

See formula III, page 25 
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(1) !/==!-. li 

(2) !/Z=l-.la;+.01a:2 

(3) y'=l-.lx+.01x^-.001x' 

(4) !/2=l-.li+.01i2-.001a;'+ 

.OOOli* 

(5) !/2=l-.la:+.01x2-.001a;'+ 

.OOOli^-.OOOOli* 

(6) v2=l-.li+.01x'i-.001i'+ 

.0001i*-.00001x*+.000001i* 

See formula IV, page 25 
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Fig. IV. 
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(1) »=(.5)^ 

(2) i,= (.6)^ 

(3) !/=(.7)^ 

(4) j,= (.S)^ 

(5) y=(..f))'' 

(6) !, = (.95)^ 

(7) y=.Q9)' 

(8) !/=(1.01)* 

(9) !/=(1.0S)* 

(10) 2,= (1.1)^ 

(11) !/=(1.2)» 

See formula, V, page 27 



Fig. V. 
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(1) !/=2-(.5)« 

(2) j/=2-(.6)* 

(3) !/=2-(.7)* 

(4) 2/=2-(.8)» 

(5) 2/= 2- (.85)=' » 

(6) !/=:2-(.9)'i'" 

(7) !/=2-(.95)* 

(8) !/ = 2-(.97)* 

(9) j/=2-(.99)* 

(10) !/=2-(1.01)* 

(11) y=2-(1.03)* 

(12) !/ = 2-(1.05)* 
(13)j/=2-(1.07)* 
(14) j/=2-(1.08)* 

See formula VI, page 28 



Fig. VI. 
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(1) Iogj(=l-.5(.5)^ 

(2) log!/=l-.5(.6)^ 

(3) log!/=l-.5(.7)* 

(4) Iog!/=l-.5(.8)* 

(5) logH=l-.5(.9)* 

(6) log!/=l-.5(.95)=' 

(7) log!/=l-.5(.98)* 

(8) log!/=l-.5(1.02)* 

(9) logy=l-.5(l.l)* 

(10) Iog!/=l-.5(1.2)»' 

(11) Iog!/=l-.5(1.3)i» 

(12) log!/ = l-.5(1.5)* 

(13) logz/=l-.5(2)^ 
base =10 

See formula VII, page 32 



Fig. VII. 
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(1) !/=2-.01x-(.5)* 

(2) i/=2-.03i-(.S)* 

(3) »=2-.05i-(.5)* 

(4) !/=2-.08i-(.5)* 

(5) j/=2-.li-(.5)* 

(6) j/=2-.12a:-(.5)* 

(7) »=2-.15i-(.5)'' 

(8) w=2-.2i -(.5)* 

See formula VIII, page 33 
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Fig. VIII 
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(1) ],=10-81-.36x + .03a;" 

(2) j,=10.54-.2te + .02i2 

(3) 2,=lO-27-.12H-.01i2 

(4) j,=lo-135-.06a: + .005i2 

(5) j/=10--135 + .06i- 005i2 

(6) j,=10--54+.24i-.02i! 

See formula IX, page 37 
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Fig. IX. 
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(1) K=(1.01)^ (1.05)'l-2)^ 

(2) ^=(1.01)* (1,06)(1-1®^ 

(3) j/=(1.01)== (1.05)(1-15)^ 

(4) a=(.5)^ (2)(l-2«'^ 

(5) a=(.5)^ (2)<1.23)^ 4:r 
.(6) !/=(.5)^ (2) (1-2)^ 

See formula X, page 37 



Fig. X. 
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10 U 18 a: 



(1) » = 



(2) y= 



(3) y 



(4) !/= 



(5) y 



.2-.lx+.05x2 

X 

.2-.li+.07l2 

X 

.2-.lx+.\x^ 

X 

.2-.ll+.2x2 

X 

.2-.lx+Ax^ 
See formula XI, page 3S 



Fig. XI. 
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Fig. XII. 



(1) »=5i-* ' 

(2) j/=5i-2, 

(3) y=5x-^ 

(4) !/=5l--l 

(5) u=5x--^ 

(6) y=5x--* 

See formula XII, page 42 
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(1) i/=l+Ioga;+.l log^i; 
^=—1.5 (min.) when log a:=~5 

(2) !,= H-log i+.Ollog' a:; 

2/=— 24 (min.) when log x=— 50 

(3) 2/=l+.21ogi+.31og'i 

(4) y=l— log i+log' I 

(5) !/=!— log a;+.51og^a; 

See formula XIII, page 44 
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aX-Ur-l+.OOSx'-'' 

(2) j,=l+.007il-6 

(3) !/=H-.006xl-5 

(4) 2/=l-.002a:2 

(5) !/=l-.003i2.1 

(6) !/=l-.004i2-Z 

See formula XIV, page 45 



Fig. XIV. 




(1) !/=(2,0) lO-'Ol^^ 

(2) j/=(1.6) io--<'2i*-'^ 

(3) »=(l.-2) 10 --OS^^-* 

(4) !/=(1.0) 10 --O*^'-^® 

(5) 2/= (0.8) io--05il-2* 

(6) !/=(0.6) lO-'OS^'-^^ 

See formula XV, page 49 



Fig. XV. 
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(1) (a-2) (i-.5) = -l 

(2) (J/-2) (s+.75) = -1.5 

(3) (!/-2) (3:+l) = -2 

(4) (»+.l) (i+4)=8.1f 

(5) (I/+.1) (1+3) = 6.3 

(6) (B+.l) (a:+2)=4.2 

See formula XVI, page 53 
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Fig. XVI. 
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(1) y=l 101+24 

-12 

(2) j,=i 10a:-24 

2 

(3) !/ = J 101+2 

2^ 

(4) !/=Ti) lOa: 

.5 

(5) y lOi+l 

See fcirmula XVIo.page 56 



S I 8 

Fig. XVId. 



APPENDIX 



141 



y. 


















































^ 


-^ 




s 

1.8 


















i^ 


^ 


















--< 






ja. 


— 


■ 










,^ 


^ 


^ 








(3) 










1.8 










































\ 




s 


\ 














^ 










i 
as 

0.G 






\ 


&, 














^ 


-^ 








^rej 


\ 


\ 






















-2C 








■^ 
























=:$ 






^ 










0.4 






















■~~ 


^ 


^ 







tO.2 
-(U 


— - 


^ 


i;::;:::^ 




^-^ 




— ' 


-" 


'' 










/■ 


■>V) 


^ 


;^ 




















/ 


^) 


/ 






















/ 


/ 
























■ 


/ 




























1 



























(1) a=.5e-<'l^+e.05i 

(2) B=2e-('5i_.5c.l2: 

(3) j/=2.25e-05j_.75e.lx 

(4) j,= 1.8e-01i_.3e-ls 

(5) y=l.Q2e--^'!-A2e--0lx 

(6) 2,=2e--<'5i-e-.01i 

(7) j,=4.2e- -21-3.56 --25^ 

(8) y=4.5e--2i-4.1e--25a; 

(9) 2/=.25e--('li-.)Ce--15i 

(10) !/=e- l*-l,le--2a: 

(11) y=.27e--'>lx-.77e--2^ 

(12) j,=e-.l2-2e--25« 

See formula XVII, page £ 



Fig. XVII. 
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(1) y=e'01aJ(1.5 cos .Ix— .5 ain .Ix) 

(2) 3/=e~-22;(1.5cos.5x— .Ssin.Si) 

(3) !/=e--l=^(.6cos .Ix+.Ssin .li) 

(4) !/=e-l^(.2 cos .31-. 1 sin .3a;) 

(5) !/=e-02i(.4cos.l6a;+.17sin.l6i) 

(6) y=.5c~-^^ sin X 

See formula XVIII, page 61 
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(1) !/=2a;-l-a;-2 

(2) !/=3i-5-2.2x-6 

(3) j/=2.3a;-8-2a;-85 

(4) !/=.li-l+.5a;-2 

(5) !/=.33x-.0012x3 

(6) j/=.2Bi-5+.05a;-8 

See formula XIX, page 65 



Fig. XIX. 
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(1) !/=15il-5(.4)* 

(2) 2/=3a;2(,S)* 

(3) »=3a:- 2(1.5)* 

(4) ^=-.53:l-5(.75)* 

See formula XlXa, page 72 
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Fig. XX. 



(1) B= 166.25-14.5 cos 1-2.75 cos 2l-10 sin i 

(2) »=167.83— 20 cos 1—4.33 cos 2^+5.5 cos 3a:-13.28 sin 1-17.32 sin 2x 

(3) j/=167.62-17.5 cos x-Z.lh cos 2i+3 cos 3i-1.38 cos 4i-12.42 sin i-lB sin 2a:- 

2.42 sin Zx 

(4) j/=167.08-17.22 cos a;-3.5 cos 2i+5.5 cos 3i-0.83 cos 4i-2.78 cos 5x+ 

0.75 cos 61-12.14 sin 1-19.05 sin 2i-sin 3i-1.73 sin 4z+1.14 sin 5x 



See formula XX, page 74 
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Approximation, Accuracy of , 114, 118, 
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Graphical Determination of, 121 
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Changes in, 131 
Curves, 10, 11, 114, 121, 122, 131, 132 

Damped Vibration, 65 
Derivative, 95 

First, 108 

Second, 108 
Differences, 11, 12, 16, 33, loi, 115, 
116 

First Order, 12, 18 

Higher Order, 12, 16, 18, 19, 22, loi 
Differential Equation, 68, 71, in 
Differentiation of Tabulated Functions, 

100, 108 

Elastic Limit, 14 

Empirical Formulas, 9, 11, 68, 90, 100 
Errors of Observation, 11, 68 
Fourier Series, 74 

6-ordinate Scheme, 74, 76 
8-ordinate Scheme, 78 

lo-ordinate Scheme, 80 

i2-ordinate Scheme, 81, 86 

1 6-ordinate Scheme, 82 

20-ordinate Scheme, 84 

24-ordinate Scheme, 85 



Gauss, 107 

Geometrical series, 11, is, 29, 42, 45, 
49, 6s 

Indicator Diagram, 114 

Integral, Definite, 114, 118, 121, 125 

Integration, Numerical, 114, 115 

3 Ordinates, ri6 

4 Ordinates, 117 

5 Ordinates, 117 
7 Ordinates, 117 

Interpolation, 11, 100, loi, 112, 114, 115 
Trigonometric, 107 

Lagrange's Theorem, 106, in 

Maximum, 113 

Mechanical Quadrature, 114 

Method of Least Squares, 11, 19, 00, 91, 

94 
Minimum, 113 
Moment of Inertia, n, 128, 129, 130 

Normal Equations, 92, 96, 99 

Observation Equations, gi, 96, 97, 98 

Ratio, n, 46 
Residuals, Squares of, 94 
Period, 74, 88, in 
Less than 27r, 74, 88 

Physical Basis, 9 

Physical Law, 9 

Planimeter, 114 

Pressure, Upward, of Water, 49 

Simpson's Rule, 119, 121, 123, 124 

Taylor's Theorem, 94, 97 

Volume, II, 123, I2S, 127 

Weddle's Rule, 119 
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